Hadron properties in the nuclear medium 
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The QCD vacuum shows the dynamical breaking of chiral symmetry. In the hot/dense QCD 
medium, the chiral order parameter such as {qq) is expected to change as function of temperature 
T and density p of the medium, and its experimental detection is one of the main challenges in 
modern hadron physics. In this article, we discuss theoretical expectations for the in-medium 
hadron spectra associated with partial restoration of chiral symmetry and the current status of 
experiments with an emphasis on the measurements of properties of mesons produced in near- 
ground-state nuclei. 
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L INTRODUCTION 

Quantum chromo dynamics (QCD), w hich is the colo r 
SU(3) gauge theory of quarks and gluons (jNambulJ 19661 ). 
is now estabhshed as the fundamental theory of strong 
interactions. The Lagrangian density of QCD reads 

q 

+ {Qi^^Qk + Qk^Qi.) (I-l) 
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FIG. 1 The running coupling constant determined from r 
decay, T decay, deep inelastic scattering, e^e ~ annihilation, 
and the Z-boson resonance shape and width (|Amsler et a/.L 
l2008l ). 
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FIG. 2 The masses of d, and s quarks at the scale k — 2 
GeV constrained by various hadro n masses using QCD sum 
rules and lattice QCD simulations (|Amsler et all l2008l ) . 



we will focus on the dynamical breaking of chiral sym- 
metry (DECS) realized in the QCD vacuum and in the 
hot-dense QCD medium by using in- medium hadrons as 
useful probes of QCD matter. 



where we focus on three light flavors q = {u, s) with 
the mass matrix m = diag(m^, m^^, m^) throughout this 
article. The quark field q belongs to the triplet repre- 
sentation of the color gauge group SU(3)c- The right 
(left) handed quark q^ = |(1 + ^^)q {q^ = |(1 - ^5)q) 
is the eigenstate of the chirality operator 75 with the 
eigenvalue +1(— 1). The covariant derivative is defined 
as Dn 



A. QCD symmetries 

Similar to the running coupling constant as(/^), the 
quark masses receive quantum corrections and become 
scale dependent, m{n). As seen from FiglSJ the current 
determination of the u and d quark masses show that 
they are about 50 to 100 smaller than the QCD intrin- 



dn + iqt^A^, with q being the strong coupling • 1 a u-i -h. 1 • ^ 

a u • oijfo\ A Aa V. • scale A^^^j^ , while the s quark mass is comparable to 



constant, being the SU(3)c generator and being 
the color-octet gluon field. The field strength tensor of 
the gluon is defined as G^^ = d^A^^-d^A^ —gfapjA^A^ 
with faf3^ being the structure constant of SU(3)c. The 
QCD Lagrangian Eq. p.ip is exactly invariant under the 
local SU(3)c gauge transformation of quarks and gluons. 

The running coupling constant g{hz) is defined as an 
effective coupling strength among quarks and gluons at 
the energy scale tz. Due to the asymptoti c free nature 
of QCD , q(K,) becomes small as k, increases ( Grossl . 120051 : 
IPolitzerl . ,2005: Wilczek, 2005). This is explicitly seen in 
the two-loop perturbation theory as 



47r/?oln(^VA2 J 



ln(ln(^VA^^J) 
ln(/.VA|oJ 



^^,/9o = (ll-iA^f)/(47r)2,/3i = (102 
Nf is the number of fiavors and A^^j^ 



IS 



where as{i^) ■■ 
f7Vf)/(47r)4^ 

called the QCD scale parameter to be determined from 
experiment. 

FiglH and Eq. p.2p indicate that the running cou- 
pling constant increases and becomes strong at low 



energies k, ^ ^qcd ^ ^00 MeV. This is the 
typical energy scale where various non-perturbative 
effects such a s the confinement of quarks and 
gluons (|Wilsonl , l2005h an d the dynamical break- 
ing o f chiral symmetry (iNambu and Jona-Lasiuiol 
ll961al lbl) (|Hatsuda and KunihirolJ 19941 ). Both effects are 
responsible for the formation of composite hadrons and 
nuclei and for the origin of their masses. In this article. 



QCD 5 

^QCD- Therefore, it is legitimate to treat itIu/Aqcd 
rrid/A^^^ as small expansion parameters. 

In the limiting case where mu,d = 0, which is called the 
SU(2) chiral limit, the QCD Lagrangian Eq. p.ip acquires 
an exact global symmetry called chiral symmetry under 
independent SU(2) rotations of the left handed and right 
handed quarks: q^ Ui,q^ and q^ ^R<?r with /7l,r 
being the global SU(2) matrices. Thus we have the exact 
QCD symmetry for rriu^ = 0, 

g = SU(3)c SU(2)l SU(2)j^ U(1)b, (1.3) 

where U(1)b corresponds to the baryon number symme- 
try corresponding to the global phase rotation, q^^^^ 
L(j^) . Although a similar phase rotation, q^ e^^q^ 



(1.2) and q^ ^ e '^^q^^ looks like a symmetry of Eq. p.ip . it 



is broken explicitly by a quantum effect known as the 
axial anomaly. Currents associated with these "symme- 
tries" are defined as = qj^t^q (the triplet vector 
current), = qj^^J^t^q (triplet axial- vector current), 
= Qlf^t^Q (baryon current), = q^^^^t^q (singlet 
axial current), where t^^i'^^s = r^/2 with being the 
Pauli matrices and = r^/2 = 1/2. The divergences of 
these currents are 



QHya 



iq[m,t'']q, (1.4) 

iq{m,t''}'j5q, (1-5) 

0, (1.6) 

ig-m75</-2gGrG^., (1.7) 



with G^^ = ^e^i^xpG^^ being the dual field strength of 
the gluon. [ , ] and { , } are the commutator and the 
anti-commutator in the flavor-space, respectively. For 
later convenience, we define the scalar and pseudo-scalar 
density as 



2.0 r 



pO 



(1.8) 
(1.9) 



From the time component of the currents, genera- 
tors of the chiral transformation are defined as Q^it) = 
/ V^{t,x)d^x and Ql{t) = J A^{t,x)d^x. Then the bi- 
linear quark operators defined above obey the following 
relations under the axial transformation Q'^ {a = 1, 2, 3); 



[Qg(t),T/;(t,x)] = +ie„fce^^(t,x), 
[QUt),Al{t,x)]=+ieabcV^{t,x), 
[Qt{t),S°{t,x)]=+iP-{t,x), 
[QUt),P°{t,x)] = -iS''{t,x), 
[QUt),S''{t,x)] = +iSabP°{t,x), 
[Qtit),P\t,x)] = -i5abS°it,x). 



(1.10) 
(1.11) 
(1.12) 
(1.13) 
(1.14) 
(1.15) 



In the past few years, remarkable progress was made 
in calculating the hadron spectra on the basis of lattice 
QCD simulations with dynamical s quarks. This 
progress was achieved partly because the supercomputer 
speed is doubled every 1.2 years and partly because of 
new simulation algorithms: The lattice QCD simulations 
for quark masses very close to the physical p oint are now 
possible in the \yilson fermion formalism ( Aoki et'oll , 
l2QQ8l : lDiirr et a/.Ll2QQ8h . Shown in FigElis an example of 
the lattice results for meson and baryon masses extrapo- 
lated to the physical quark masses using the simulation 
data in the interval, ^{mu-\-md){i^ = 2GeV) = 3.5 MeV— 
67 MeV (corresponding to = 156 MeV - 702 MeV). 
The experimental data are reproduced with 3% accuracy. 
The simulations right at the physical quark masses will 
be performed in the very near future. 



B. Dynamical breaking of chiral symmetry in the vacuum 

Even if QCD in the SU(2) chiral limit has the sym- 
metries of Eq. p.3p . the ground state of the system may 
break some of the symmetries dynamically. Let us con- 
sider the QCD vacuum |0) at zero temperature and den- 
sity. Assuming that the vacuum is Lorentz invariant and 
taking into account the fact that QCD does not allow dy- 
namical breaking of parity and vector symmetry in the 
vacuum (IVafa and Witte 3,11984'), we have the following: 
possibility of the symmetry breaking pattern; 

SU(2)l X SU(2)j, ^ SU(2)l+r ^ S\J{2)^. (1.16) 

In terms of the generators of the vector and axial- vector 
rotations, such a vacuum state is characterized as 
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FIG. 3 Light hadron spectrum obtained from lattice QCD 
simulations with dynamical s quarks in the Wilson 
fermion formalism. The hadron masses are extrapolated to 
the physical quark masses (determined by m^r, rriK and mn) 
using the data in the interval, ^{ruu + md){i^ = 2GeV) = 
3.5 MeV — 67 MeV. The spatial lattice volume V and the 



lattice spacing a are (2.9 fm)^ and 0.09 fm, respectively , 
zontal bars denote the experimental values (|Aoki et a/.l J 



Hori- 



Strictly speaking, we need to take the SU(2) chiral limit 
^n,(i after taking the thermodynamic limit V ^ 
to make the matrix elements of well defined. This is 
similar to the case of the spin system where the external 
magnetic field plays the role of mu^d- 

At this point, it is in order to mention general defini- 
tion of the order parameter. Consider a symmetry group 
Q and its generator Q. If there is an operator ^ such that 
{[iQ,^])o{= (0|[i(5, <I>]|0)) 7^ 0, this expectation value is 
called the order parameter. If the vacuum is symmet- 
ric under Q, the order parameter becomes zero. On the 
other hand, if the vacuum is not symmetric under Q, 
there exists a Nambu-Goldstone boson having the same 
quantum number as <I>. Note that the order parameter 
is not unique for a given Q: one can introduce higher di- 
mensional o rder parameters in principle to characteri ze 



mensional o rder parameters m prmciple to cnaracteri ze 
the system (|Kogan et a/.IJl999l : IWatanabe et a/l,l2004l ). 

For the symmetry breaking pattern as Eq. (|I.16p , Q 
is identified as Qg, and a simplest choice of <I> is P^. 
Then, it leads to the order parameter {S^)o. Recent lat- 
tice QCD simulation of the chiral condensate using over- 
l ap Dirac ferniion wi th dynamical u^d^s quarks indicates 
(|Fukava et all l2009h 



1/- 
-(uu- 



-dd)o = -(242(04)(i{^) MeV)^(T18) 



Q«|0)=0, Qg|0)7^0. 



(1.17) 



where the renormalization scale is taken to be = 2 GeV 
with the statistical and systematic errors in parentheses. 
This result implies that the QCD vacuum is the Bose- 
Einstein condensate of quark— anti-quark pairs {qq)o = 
((^L^R +^r^l))o and has the power to change left handed 
quarks to right handed quarks and vice versa: Namely 
the condensate induces a dynamical quark mass. Since 
quarks are confined, it is not possible to isolate a sin- 
gle quark to measure the dynamical quark mass. Nev- 
ertheless, there is indirect evidence that the quarks in- 
side hadrons have an effective mass (constituent quark 
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FIG. 4 Effective quark mass M as a function of the Euclidean 
momentum q obtained from lattice QCD simulations with- 
out dynamical quarks (quenched) and with three dynamical 
quarks (unquenched) ([Bowman et a/.L l2005l l . 



mass) M ~ 350 MeV from the phenomenological de- 
scription of hadrons. The effective quark mass near zero 
Euclidean momentum in lattice QCD simulations with 
Landau gauge fixing leads to a similar value as shown in 
FigH 



Kl 1 
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FIG. 5 Light scalar (5'), pseudo-scalar (P), vector {V) and 
axial- vector [A) mesons. Their spin and parity are denoted as 
J^. Threshold of the tttt, and KK decays are also shown 
(|Yagi et ali l2Q05l ). 



partner at the same mass, which are the direct evidences 
of DBCS. Such non-degeneracy is also seen in other chan- 
nels, e.g. cj, K* and (j) as illustrated in Fig El 



C. Chiral symmetry and hadron spectra 

The Nambu-Goldstone (NG) bosons associated with 
the DBGS of SU(2) chiral symmetry are nothing but 
the charged and neutral pions. Moreover, the par- 
tially conserved axial current (PC AC) relation, Eq.(|L5]), 
l eads to the Gell-Mann — Oakes— Renner (GOR) relation 
( Gell-Mann et a/.LIl968[ ) which relates the pion masses to 
the quark masses as 



c2 2 
2 2 



-m(uu + dd)o + O(m^), 
-(muUU + mddd)o + 0{m^). 



(1.19) 
(1.20) 



Here m={m^ + ^d)/2 is the averaged mass of u and 
d quarks, fj^^ = 92 A MeV) is the pion decay constant, 
and m^± 140 MeV (m^o ^ 135 MeV) is the charged 
(neutral) pion mass. Using these values together with 
the quark masses in Fig|2l we obtain the finite chiral 
condensate comparable to Eq. p.lSp . 

Further experimental evidence of DBCS is obtained 
from the observed meson spectra. If chiral symmetry 
is not broken, QslO) = 0, vacuum expectation values of 
all the the commutators, ([Qs", • • • [QsS ^]] • • •])o, 
should vanish for an arbitrary operator <l>. The con- 
traposition of this statement with n = 2 and ^ = 
6'^(x)5'^(^), V/^{x)V^{y) leads to a statement that DBCS 
must occur if the correlation functions of the chiral part- 
ners are not degenerate. Namely, 



(5»(x)5°(y) - 



-A;ix)Aiiy))oy^O 



Qt\0) ^ 0, (1.21) 
.gg|0)7^0. (1.22) 



Experimentally, the pion (the pseudo-scalar meson) does 
not have a scalar partner at the same mass, and the p- 
meson (the vector meson) does not have an axial-vector 



II. CHIRAL SYMMETRY AND IN-MEDIUM HADRON 
SPECTRA 

Connections between properties of the QCD vacuum 
and hadronic correlation functions as discussed in Sec lLBI 
and Sec jLCI can be generalized to QCD at finite temper- 
ature and density. In this section, we will summarize 
such theoretical connections with special emphasis on the 
pion, the scalar meson, and vector mesons in the medium. 



A. Chiral condensate in the medium 

Let us now consider how the simplest chiral order pa- 
rameter (qq) changes its value inside the hot and/or dense 
medium. Exact formula for the in-medium chiral conden- 
sate in terms of the QCD partition function Z at finite 
temperature T and the baryon chemical potential fi is 
given by 



where 

Z(T,m) 



Tr 



qqe 



-Kqcu/T 



-Kqcjd/T 



drua 



^P{T,fi)V/T 



(IL23) 



(IL24) 



i^QCD = H^cD + J {Qrnq - q^fiq)d^x, (IL 



25) 



with Hqq^ being the QCD Hamiltonian without the 
quark mass term. 
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FIG. 6 Normalized chiral condensate A = [{uu)t — 
{m/rris){ss)T]/[{uu)o — (?fi/ms)(ss)o] as a function of T for 
two different lattice spacings, a = 0.24 fm (A^r = 4) and 
a = 0.17 fm {N^ = 6), calculated by the lattice QCD simula- 
tions with dynamical u. d, s q uarks in the staggered fermion 



tions with dynamical u. a, s q u; 
formalism (|Cheng et aa 120081 ). 



1. Finite temperature 



Eq. pL23p can be evaluated analytically in some spe- 
cial cases. For example, in the two-flavor system (A/'f = 2, 
i.e., rriu = rrid < oo, = oo) at low T with /i = 0, the- 
mal pions are the dominant contributor to the pressure. 
In the leading order of t he virial expansion by t he pion 
number density, we have (|Gerber and Leutwvlerl , [l989) 



{qq)o 



Nf^2 



6/r 



e(T) = B^imjT), 



(11.26) 
(11.27) 



^Bi(mT,/T) with ub be- 



ing the Bose-Einstein distribution and Sk = y k 
Eq. pi.27p shows a clear tendency that the magnitude of 
the chiral condensate decreases as T increases. However, 
for T > 150 MeV, the interaction among pions and the 
contribution from other mesons become important and 
the estimate based on Eq. ()II.27p is not reliable. The two- 
flavor system at extremely high T with /i = can be also 
evaluated because the thermal quarks and gluons are the 
dominant contributor to the pressure due to asymptotic 
freedom. In the leading order of as(/^ ~ T), we have 
{qq)T — \mqT'^ {q = u^d), so that the chiral condensate 
vanishes in the chiral limit. Lattice QCD simulations at 
finite T in Fig l6] indeed show a decrease of the normalized 
chiral condensate for u, d quarks w i th a r apid crossover 
around T 200 MeV (|Cheng et a/.Ll2008h . 



Pq = and fiq = -^s with e being the energy den- 

sity of the system: 



{qq)p {qq)o 



d 

dmrin 



(11.28) 



where p{= ^qpq)^ A and E are the total baryon 
density, the total baryon number, and the total en- 
ergy, respectively. There is an alternative derivation 
of this formula usi n^ the Hellmann-Feynman theorem 
(ICohen et a/.LIl992h . 

For two-flavor nuclear matter with equal num- 
bers of protons and neutrons, the leading order of 
the vir ial expansion in t e rms of the baryon d ensity 
reads ([Coh en et al). 'l992l: iDrukarev and Levinl . Il99ll : 
iHatsuda an d Lee, 1992) 



{uu + dd)p 
{uu + dd)o 

{ss)o 



1 



y~P? — 9P Fi[kF/mN)' 



7>r 



pFiikF/ruN), (n.29) 
(11.30) 



Here ctttAt = m{N\uu 
and y = 2{ss)n/ {uu 
tent of the nucleon. 



\- dd\N) is the ttN sigma-term, 
\- dd)N is the strangeness con- 
Also, /72™^(fcF - |fc|) 



Q^Fi(kY /rriN) with /cp being the Fermi momentum and 
Ek = \lk^ + m- 



. The low density expansion of Fi reads 

'^2,9 



56 



• (jHatsuda et all Il995h . 



Using the empirical values, cr ^/y = 45 =b 10 MeV and 
7/ = 0.12 - 0.22 fCasser etM^ Il99ll; iHatsuda and Led . 
Il992h the right hand side of Eq. (|IL291) gives almost 
a 35 % reduction of {qq)p at nuclear matter density 
po = 0.17 fm~^. This leading order result together with a 
calculation of higher order corrections on the basis of in- 
medium chiral perturbation theory are s hown in Fig 13 for 
nucle ar matter and neutron matter (Kais er and Weisa 
I2009D . 

The change of the chiral condensate induced by the 
strong electric and magnetic fields and strong color elec- 
tric and magnetic fields is also an interesting subject 
which may be relevant to the physics of relativistic heavy 
ion collisions and th e structure of conipact star s with 
high m agnetic field (Klevansky!, Il992l : iMiranskvl . l2002l : 
ISuganu ma and Tatsumi, 1991). 

The experimental detection of the change of chiral con- 
densate in the hot and/or dense medium is one of the 
most interesting challenges in modern hadron physics. 
Possible tools are hadron-nucleus and photon-nucleus re- 
actions, heavy-ion collisions, and deeply bound mesic 
atoms and mesic nuclei, which are summarized in the 
later chapters of this article. 



2. Finite baryon density 

A useful formula for the chiral condensate at fixed 
baryon density p with T = is obtained from Eq. pI.23P 
by using the thermodynamic relations; P-\-e = PqPq^ 



B. Spectral functions in the medium 

In the hot and/or dense environment, all the hadrons 
including the pion undergo spectral changes due to their 



6 



m =135MeV 



Isospin-symmetric nuclear matter 



correlation, 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 



0.8 
2 0.6 

1 0.4 

o 

0.2 
0. 



' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' 


' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 . 




^^^^^m^=135MeV ■ 












linear ^ ^^^s. \ 


Neutron matter 

" , , , , 1 , , , , 1 , , , , 1 , , , 


. \ .... \ .... \ \' 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 

FIG. 7 Comparison of the linear density approximation of the 
chiral condensate (dashed lines) and the results with higher 
order terms originating from nucleon-nucleon correlations cal- 
culate d by the in-medium chi ral perturbation theory (solid 
lines) ([Kaiser and Weisd . l2QQ9l ) . 



strong interactions with the medium. Therefore, it is 
not enough to talk about the "mass" and "width" of the 
hadrons, but we need to study the hadronic spectral func- 
tions. For mesic resonances coupled to the composite 
operator 0{t^x) = q(t^x)Tq(t^x) with F being an arbi- 
trary combination of Dirac and flavor matrices, we have 
the spectral decomposition of the retarded correlation as 

G^{uj,k) = i¥.T.{R{0(t,x)0\{)))), (1131) 

— {uo ^ iS)' 



-du' 



(11.32) 



where F.T. stands for the fourier transform, 

R{A{t,x)B{t',y)) = e{t - f)[A{t,x),B{t\y)] is 
the retarded product of the operators A and 5, and (•) 
implies the expectation value at finite T and /i. The 
spectral function, p^iuo^k)^ has all the information of 
the states having the same quantum numbers with the 
operator O. In particular, and V^^ are the relevant 
operators for studying the in-medium pion and the 
p-meson, respectively. 



1. In-medium pion 

The Nambu-Goldstone theorem which guarantees 
massless pions in the QCD vacuum in the chiral limit 
holds also in the medium. Indeed, by considering the 



n = J d^xd^{R[A^^{t,x)P\0)]), (n.33) 
the following sum rule can be derived (lYagi e.t all \m)^: 

L 



2m 



(11.34) 



where we considered a two-flavor system (m = = 
rrid < oo,ms = oo) for simplicity, p^^ corresponds 
to the spectral function associated with n^^(x) = 
(R(P^(t, ic)P^(O))). To have a non-zero chiral conden- 

0, 



sate in the r.h.s. of Eq. (|II.34p in the chiral limit m 
the spectral function must have a pole, p^^(cc;, 0)|m^o ^ 
S^^[CS{ij'^ — am)-\-- • •], so that m in the numerator is can- 
celed by this pole in the chiral limit. This means nearly 
massless pions exist even in the medium as long as DBCS 
takes place. 

In the leading order of the virial expansion at T 7^ 
with /i = 0, the self-energy of the in- medium pion 
is dictated by the forward pion-pion scattering ampli- 
tude which vanishes in the chiral limit. In this case, 
the pion is still a real pole with the mass m^(T) and 
the decay constant fj (T) fOoitv and Leutwvle E Il989l : 
iPisarski and Tvtgatl . Tl996l : pToublam 19971 ): 



u J 



i-e(T). 



(11.35) 
(11.36) 



Here, the in- medium pion decay constant fl{T) is defined 
by the residue of the pion pole of the correlation function 
of A^^Q. Since the Lorentz symmetry does not hold in 
the rest frame of the medium, a difference arises between 
the temporal residue fl and the spatial residue in 
general, although they are equal in the leading order of 
the virial expansion at finite T. Combining Eq. pI.36P 
and {qq)T in Eq. (|II.26p . the GOR relation turns out to 
hold i n the dilute pio n gas at low T ( Pisarski and Tvtgatl . 
[i996[| Toublaij ,E993): 



{fi{T)m^(T)y 
m{uu + dd)T 



(11.37) 



which was original ly noticed in the Nambu— Jona- Lasinio 
model at finite T (|Hatsuda and Kunihirol . fl987 bl). 

In the case p 7^ with T = 0, the in- medium pion 
properties at low baryon density is determined by the 
pion-nucleon forward scattering amplitude whic h is dic- 
tated by several low energy constants (Meissn er et alV 
l2QQ2l : [rh orsson and Wirzbal , Il995[ ). For the symmetric 
two-flavor nuclear matter, we have 



m^{p) 



1 + ^1 2ci 



C2 - C3 



Stun, 



p(II.38) 
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= ^ + j2{-- + -^-^)P^ ("-39) 



Using Eq. pL29p and a relation cJt^n — — 4cim^, the GOR 
relation is shown to hold at low density: 



rh{uu + dd)p 



-1. 



(11.41) 



With the empirical values, ci = -0.81 ± 0.12 GeV"\ 
C2 = 3.2 ± 0.25 GeV"\ cs = -4.70 ± 1.16 GeV"\ 
= 92.4(3) MeV and qa = 1.2695(29) MeV, one finds 
fi{p)/U = l-(0.26±0.04)(p/po), mp)/U = 1-(1.23± 
0.07)(p/po), {qq)p/{qq)o = 1 - (0.35 ± 0.09)(p/po). Note 
that, in asymmetric nuclear matter, there is splitting be- 
tween m^- ^ and m^o: For example, N/Z = 1.5 at 
P = Po = 0.17fm~^, there is approximately a +18 MeV 
shift fo r TT", a -12 MeV shift for 7r+, and a +2 MeV shift 
for TT^ (|Meissner et a/.L [20021 ). A possible problem of the 
linear density formula Eq. pi.40p is that vanishes even 
below nuclear matter density. 

In the leading order of baryon density and in the SU(2) 
chiral limit, one can formulate two different representa- 
tions of the in-medium change of the chiral condensate 
in terms of the physical observables: 



{qq)p 
{qq)o 



zy\p) 



hi 



h{p) 



1/2 



(11.42) 
(11.43) 



The first one is the in-medium GOR relation 
Eq. pl.4ip . while the second one is a combina- 
tion of the in-medium Tomozawa- Weinberg relation 



(jKolomeitsev et l2003bl ). bi/bi{p) ( fi(p )/U?. and 
the in-medium Glashow- Weinberg relation (jjido et all 

Uli), {qq)p/{qq)o ^ zV\p) x Here hr{p) 

{hi) is the isovector pion- nucleus (pion-nucleon) scatter- 
ing length in the chiral limit and Z^{p) is the in- medium 
change of the pion pole residue of the correlation func- 
tion of the pseudo-scalar operator P^{x). The slope of 
Zt^{p) as a function of p at low density is related to the 
isoscalar pion-nucleon scattering amplitude, while hi{p) 
is related to the energy levels of the deeply bound 7r~- 
atom. Experimental data at p < po indicate Z-j^{p) < 1 
and bi/bi{p) < 1, so that the chiral condensate indeed 
decreases at finite baryon density. 



2. In-medium scalar meson 

The light scalar-isoscalar meson has been custom- 
arily called the a. Since it has the same quantum 
number as the vacuum, a is analogous to the Higgs 
boson H in the electro-weak (EW) theory. The a 



may be interpreted as the excitation associated with 
the amplitude fluctuation of the chiral c o ndensa te 
aa) (iNambu , 196QI: Nambu and Jona-Lasinio, Il961allb ) 



Gell-Mann an d Lev 
1982; WeinberJIm 



1960l ) ( Delbourgo and ScadronI , 
However, there is a marked dif- 
ference between H and a: In the EW theory, the NG 
bosons associated with spontaneous symmetry breaking 
SU(2) X U(1)y U(l)em are absorbed into the gauge 
bosons, while the NG bosons in QCD (the phase fluctu- 
ation of (qq)) are nothing but physical pions. Therefore, 
<j is allowed to have s-wave decay into two pions as long 
as rria- > 2m^. Thus a should be a very broad resonance 
even if it exists. 

Because of the above reason, it has been long debated 
whether there is unambiguous experimental evidence of 
such a light and broad resonance in tt- tt sca ttering, 7 — 7 
collis ion, heavy meson decays, and so on (|Penningtonl . 
|2007[ ). Recently, an analysis based on the model inde- 
pendent Roy equation for the partial wave amplitude in 
the scalar-isoscalar channel, t5=o('^)' been carried out 
using precise inputs of the tt — tt scattering lengths ob- 
tained from chiral perturbation theory. The mass and the 
width of a corresponding to the s econd sheet pole are 
then deduced w ith high accuracy ( Caprini et all l2006l : 
lLeutwvlerlJ2008h : 

= Ultl^ MeV, = 544^^^ MeV. (11.44) 

Although the existence of a is established, its quark- 
gluon structure is still unkn own and is actively studied 
theoretically, experimentally ( Penningtonl 2007) and also 
numerically in lattice QCD simulations (|Kunihiro et all 
l2009l : lPrelo^^l2008D . 

The medium modification of a has not been es- 
tablished yet even at low temperature and den- 
sity unlike the case of the pion. Nevertheless, 
we may expect from general grounds that there 
would be a partial degeneracy between a and tt 
if the system approac hes to the point of chiral 
symm etry restoration (Hatsuda and Kunihiro . Il985 



1987albl) (iBernard alj [1987D (1 Chiku and Hatsudal. 
1998br [Hatsuda and Kunihirol . I2OOII : iHatsuda et al 
19991). At finite T, such a chiral degeneracy can be de- 



tected, e.g., by the thermal hadronic susceptibilities asso- 
ciated with the operator 0(r, x) defined in the Euclidean 
time r; 



Xc 



J^^ J d^x{O{r,^)O'^{0,0))T, (11.45) 



dco' 



iPo (^) 



(H.46) 



Shown in Figl8]is the lattice QCD simulation of Wl/Xo 
with dynamical quarks in two-flavor (Karsch, 200j|). One 
can see the degeneracy between the susceptibilities in the 
a channel and the tt channel as T increases to the left. 
Also, ^J\|Xo ^ is smaller than that for ao (scalar- 
isovector meson, traditionally called S) at low T, which 
indicates that the spectral strength in the a channel has 
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FIG. 8 Thermal susceptibilities in three different channels 
(7r,cr, and ao) for two-flavor QCD with staggered fermion on 
the 8^ X 4 lattice with mu,dCi = 0.02. The vertical (horizontal) 
axis denotes y/l/Xo (^^^ lattice coupling : (|Karschl , 
I2OO2I ). Low (high) T corresponds to the left (right) of the 
figure. 



more weight in the low frequency region (a is lighter than 
ao) as can be seen from Eq.(|lL46]). Splitting between ao 
and TT even when a-n degeneracy is realized at high T 
reflects the explicit breaking of U{1)a symmetry. We 
note that the close relevance of the scalar-isoscalar sus- 
ceptibility at finite baryon densi ty to the nuclear mat- 
ter p roperties is also pointed out ( Ericson and Chanfravl , 
I2QQ7D . 



3. In-medium vector meson 



Unlike the case of Eq. ()IL34p which relates the 
chiral condensate {qq) and the spectral function in 
the pion channel, no such relation is known in the 
vector channel. Still, one can derive useful rela- 
tions by using the in- medium generalization of the 
QCp sum rules ( Hatsuda et all 1 19931 : iHatsuda and Le 
Il992h. In particular, t h e in- m edium Weinberg relations 
([Kapusta and Shurvakl . ll994D for the vector and axial- 
vector spectral functions at zero spatial momentum read 



) 



(11.47) 
(11.48) 







^ 47r 

dwWip.iu;) - p,H) = -— a,(04,), (11.49) 

^ 



where O^q = 0^^^ 20^^ with 
4 



O 



Note that there are longitudinal and transverse spectral 
functions in the medium, but they coincide at zero spa- 
tial momentum, so that such a distinction is not made 
in the above formula. As is obvious from Eq. pi.50p . the 



DBCS in the vector and axial-vector channels is mani- 
fested as the higher dimensional four-quark operator and 
not by the simple bilinear operator qq = q^^q^ -\- q^qi^'- 
The in- mediu m changes o f {Oaq) and {qq )'^ are different 
in general (^Eletskvl . Il993l : iHatsuda et aU , [1993,). At fi- 
nite T with zero baryon density, it has been proven that 
there is no exotic phase in whi ch {qq) = and {O/^q) 7^ 
take place simultaneously (Kog an et a/.IJl999l ). However, 
such a phase is not ruled out at finite baryon density 
and is indeed realized in the color superconducting phase 
(IHatsuda et a/.Ll2008h . 



The spectral modifications of vector and axial-vector 
channels at low T with zero fi are realized as a mixing 
of the two channels due to thermal pions. In the leading 
order of the viria l expansion in the two-flavor system, one 
finds (jPev et all [T99Q) 



Pv 
Pa 



(1 
(1 



e(T))p; 
e(T))p^ 



(H.51) 
(H.52) 



with 0(T) given by Eq. pi.27p . The spectral functions in 
the vacuum are measured experimentally (see, e.g., Fig|9] 
in the vector channel). The above mixing formulas show 
that the pole positions of the correlation functions do 
not change at low T, while the pole residues are mod- 



ified as UliT)/fp) 



e(T). Note also that these 



formulas satisfy the Weinberg sum rules. In fact, the 
T-dependence of {OAq)T in the r.h.s. of Eq. pi.49p calcu- 
lated by using the soft pion the orem coincides with t hat 
obtained from EQs. pi.51III.52p (jHatsuda et a/.LIl993h . 

The density dependence of the four-quark condensate 
{PA.qi p is not known precisely. In the leading order of 
the virial expansion in terms of the baryon density, we 
have {p)p ^ ((9)o + {0)nP' The nucleon matrix el- 
ement oi O = 04q corresponds to higher twist terms 
in the deep inelastic lepton-nucleon scattering, but the 
value is still uncertain. A crude approxi mation origi- 
nally made was a factorization ansatz (Hats uda and Lea 
Il992h : {qiqjqkqi)N oc {qq)o{qq)N with appropriate Fierz 
coefficients. It is not obvious, however, whether this esti- 
mate is accurate enou gh and further studies are necessary 
(jThomas et a/.Ll2007h . 



C. Dynamical approaches to in-medium hadrons 

Although there are numerous attempts to relate 
hadronic sp ec tral fu nctions to (qq) in the medium 
(lAlam e^g/l l2QQll: ICassing and Bratkovskavai 1 19991 : 
iMosei I2QQ8I : iRapp and Wambachlj2QQQl ), no rigorous re- 
lations have been established yet except for the pion. In 
the following, we briefly outline various theoretical ap- 
proaches that have been attempted so far. 



1. Naive quark model 

Assuming that the constituent quark mass M origi- 
nates mainly from DBCS according to the idea of Nambu 



9 



resonance 




-f-T 



continuum 



sLGeV' 



FIG. 9 Spectral function in the isovector channel in the vac- 
uum with s = o;^) obtai ned from the e^e~ annihila- 
tion into even numbers of pions (jKwon et al 



and Jona-Lasinio, and assuming further that the vector 



meson mass follows the additive rule 



2M, one 



may expect a reduction of associated with the par- 
tial restoration of chiral symmetry. Such a shift could be 
detected throu gh the decay o f the neutral vector meson 
into dileptons ( Pisarskil , Il982[ ). 



on various models ( Eichstaedt et all l2QQ7l : iKwon et"a/] , 
I2QQ8I : iThomas et all \200^ . 



4. Hadron mass scaling 

It was conjectured that the masses of light vector 
mesons (p, co) scale u niversally as a function of density 
and/or temperature ([Brown and Rhd Il99l[ ). Near the 
chiral restoration point, {qq)/{qq)o ^ 1, the scaling law 
reads 



{qq)o' 



(11.53) 



where m* denote the pole mass of in-medium vector 
mesons (|Brown et a/.L 120081 ). The theoretical foundation 
of such a scaling law may be obtained by an approach 
in which the vector mesons are considered to be the 
gauge bosons associated with the hidde n local symmetry 
of the chiral effect i ve Lagran^an dHarada et all 12002 ; 



Harada and Sasaki 
20031 : iHidaka et al. 




20061 : iHarada and Yamawaki . 



2. Nambu— Jona-Lasinio model 



5. Bag model 



As a field-theoretical model to treat the meson proper- 
ties in the medium beyond the simple additive rule, the 
Nambu— Jona-Lasinio (NJL) model at finite temperature 
and density has been studied. In particular, the spectral 
degeneracy between tt and a in hot and /or dense mat- 
ter was expli ci tly demonstrat e d (iHatsuda and Kunihird 
119851 . Il987al ) (jBernard et all Il987h . Further progress 
along these lines and a similar dynamical model based on 
the Dyson- Schwinger equatio n can b e seen in ([Buballa . 
20051: IHatsuda and Kunihirol . [1994"; 'Kl evanskvl. Il992 : 



Vogl and Weisd . 1991 ) and in (^Roberts and Schmidt 
20001) . respectively. 



The bag model is a phenomenological approach in 
which quarks and gluons are confined in a "bag" in- 
side non-perturbative QCD vacuum. In the nuclear 
medium, the non-vanishing expectation values of the 
scalar-isoscalar meson (a) and the time-component of the 
vector-isoscalar meson {coq) develop, so that they act as 
effective scalar and vector potentials on the quarks in- 
side the bag. In such a quark-meson coupling model, 
in-medium baryons composed of three valence quarks 
feel both scalar and vector potentials with opposite sign, 
while the in- medium mesons composed of quark and anti- 
quark feel only the scalar potential and obey a universal 
scaling law dictated by the density dependence o f {a) 
(jSaito and Thomasl . Il995l : ISaito et all Il997bl . l2007l ). 



3. QCD sum rules 



The QCD sum rule (QSR) is a method which can 
relate the hadronic spectral functions to the QCD 
condensates through the ope rator product expans ion 
and the dispersion relation ( Shifman et all Il979al lbl). 
This approach has been generalized to attack the 
problems of hadron properties in the h ot /dense medium 
( Bochkarev and Shaposhnikovl. 1 19861 : IHatsuda et all 

and further theo- 



retical elaborations were made (Asakawa and Ko. 


1993; [Jin and Leinweber 


. 19951: iKlind et al. 


1997 


Koike and Havashisaki, 


19971: iLeuDold et al. 


1998 



gppert et all I2006D . A major difference of the in- 
medium QSR from the in- vacuum QSR is that there 
arise Lorentz-tensor condensates. The weighted average 
of the spectral function (p)^ = J duo p{u;)W (co) obtained 
from the in-medium QSR gives useful QCD constraints 



6. Hadronic models 

There are purely hadronic descriptions of the in- 
medium vector mesons. In the Walecka type models, 
the one-loop self-energy of vector mesons in nuclear mat- 
ter receives contributions from the low-energy particle- 
hole (p-h) excitations and the high energy nucleon— anti- 
nucleon {N-N) excitations. The former gives the stan- 
dard plasmon effect which increases the pole mass of 
the vector meson, while the latter tends to decrease the 
pole mass: These effects are simply understood from 
the quantum mechanical level repulsion. The net ef- 
fect with a standard parameter set of the Walecka model 
shows tha t the N-N effect wins and th e pole mass 
decre a ses (Caillon and LabarsouQud . Il993l: IJean et all 
\1994 iKurasawa and Suzukil . Il990l : ISaito et all Il989 : 
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[Sh iomi and Hatsudal . [T99l. For more details of such a 
model, see the review (iHatsuda erd] , Il996h . The ac- 
tual vector meson in nuclear matter receives not only the 
nucleon effect but also effects from resonances such as A 
and . They modify the spectral structure non-trivially 
and cause s pectral shift, spectra l ^broadening, and even 
new peaks ("Asakawa et al\ 1992; Chanfrav and S chuck , 
1993; Friman and Pirner, 1997; Herrmann et al\, 1 1993 : 



Post all l2QQ4l : iRapp et all Il997l : iRaPD and Wambach 
2QQQI ). 



7. Chiral effective theories 

Spectral change of the broad resonance a in the hot 
and/or dense medium is one of the interesting signals 
of ch i ral sym metry restoration ( Hatsuda and Kunihirol . 
119851 . Il987al ). An example is the enhancement of 
the spectral function in the sca lar- isoscalar channel 
near t he 2tt threshold at finite T ( Chiku and Hatsudal , 

Il998h a nd at finite b aryon 



Volkov et al 



density (Aouissat et all 1 19951 : iHatsuda et all 1 19991 : 



ISchuck et all\l9S^ . Since a and p are both resonances in 
the TT-TT system, it is necessary to take into account mul- 
tiple scattering of the pions inside the medium to study 
their spectral structures. For this purpose, unitarization 
of the in-medium tt-tt scattering amplitude of the chi- 
ral Lagrangian on the basis of the N/ D method and the 
inverse amplitude method has been extensively explored 
(ICabrera a/.l.l2QQ8l .'2QQ5: Fernandez -Fraile e^"aIl.l2QQ7 
Jido et all l2QQll : lOller et aL 2QQQ : lYokokawa e7 



2QQ2I ). The large expansion of the 0{N) symmet 



ric scala r mod el has als o been studied as a toy model 
(fHidaka et all [2004 : .Patkos et a/.l . l2002h . 



8. Lattice QCD 

The best quantitative method to study the medium 
modification of hadrons should be the lattice QCD 
simulations. The maximum entropy method (MEM), 
which allows us to extract the in-medium spectral func- 
tions from the Euclidean correlation functions measured 
on the lattice has been proposed and tested in the 
quen ched lattice QCD simulat i ons at finite tempera- 
ture ( Asakawa et~al\ . I2OOII I2OO3I : iKarsch and Laermannl . 
|2003[ ). Its application to the lattice data with dynam- 
ical quarks has been started at finite T although it 
i s still limit e d to heavy quarkoniums at the moment 
( Aarts et all |2007[ ). The application of this method to 
QCD at finite density is still not attainable due to the 
sign problem which invalidate s the importan c e sarn pling 
in Monte Carlo simulations (|Murova et al I I2OO3D . In 
the mean time, it would be useful to study in- medium 
hadronic properties by analytic lattice approaches such 
as the strong coupl ing expansion at finite T and /i 
(jOhnishi et a/J . l2008D . 



III. PSEUDOSCALAR MESON: tt IN NUCLEI 
A. Theoretical background 

As we have shown in the introduction, a sizable reduc- 
tion of the chiral condensate in nuclear matter of about 
30% is theoretically expected. A possible way to detect 
this reduction is to study the in-medium pion properties 
through precision spectroscopy of deeply bound pionic 
atoms and through precision measurements of the low- 
energy pion-nucleus scattering. 

The basic tool to relate the in-medium chiral conden- 
sate and the experimental data is the pion-nucleus optical 
potential t/opt(r). The starting point is the pion propa- 
gator in asymmetric nuclear matter. 



n(cj,q). (III.54) 



Here 11 denotes the in-medium pion self-energy which is 
written in the linear density approximation as 



(III.55) 



where p = Pp -\- Pn (total baryon density), 6p = 
Pp — Pn (total isospin density), and e = +1,0,-1 for 
7r~,7r^,7r+. The off-shell isoscalar (isovector) forward 
scatte ring amplitude "^^ ^ "^ fof small co and at q = 
reads (|Ericson and Weis Il988[ ) 



r+(c.) 



J TT 



From the constraints T+(a; = m7r,0) = 47r(l + 
ttIt^ /M)aT^N with the accidentally small i soscalar scatter- 
ing le ngth a^N = (0.0016 ± 0.0013)m-i ([Schroder et all 
I2OOII ). we have (3 a^^N/Tn'^. This approximation is 
valid in 5% accuracy to p. The tt^ does not have a mass 
shift in the nuclear medium in the same approximation. 
One may introduce the energy independent optical po- 



tential as uj'^ 



ml 



the dispersion relation uo'^ 



2m-j^Uopt which is a solution of 



n(cc;) = 0. As long as 



^opt ^ and the linear density approximation is valid, 
one finds 



Sp 



opt 



mipW 



(III.57) 



This implies that there is an extra repulsion (attraction) 
for 7r~ (7r+) from the medium associated with the reduc- 
tion of the pion decay constant (Weise, 2000, 2001a). 

In the local density approximation, pp^ri Pp,n{^)i 
the Klein-Gordon equation corresponding to Eq. (|III.54p 
reads [u'^ + - - n(cj, V; pp,n(r))] ^{y) = 0, or 
equivalently near the mass shell, 

[uo^ + V' - - 2m,/7opt(Pp,n(r))] $(r) = 0,(111.58) 

with = uo — Fcoui(r) where Vcoui is the Coulomb po- 
tential between tt^ and the nucleus. There are several 
important contributions to Uo^t other than Eq. pil.57p 
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in studying the experimental data of the deeply bound 
pionic atoms and the low energy 7r-nucleus scatter- 
ing: (i) higher order terms in uj and m^r t o the 
s-wave part of and T~ (|DQring and Osetl . l2QQ8l : 
iKolomeitsev et all . l2QQ3al lbl). (ii) the two-nucleon ab- 
sorption of the pion, and ( iii) the p-wave contribu- 
tion with spatial derivatives (lErics on and Wei"s3 . 1 19881 : 
[Friedman and Gal, 2007). Their explicit forms are given 
in the next section. 



B. Pion-nucleus optical potential 



1. s-wave and p-wave parts 



Th e pion-nucleus potential (|Ericson and Ericsonl . 
1 19661 ) is composed of the s-wave and p-wave parts: 



f/opt(r) = U,{r) + [/p(r), 
Us{r) 



27r 

[6(r)+e2Bop'(r)], 

27r - 



(III.59) 
(III.60) 



f^P W = V • [c(r) + £2"'Cop2(r)]L(r)V, (III.61) 



m 



with 



b{r) 
c{r) 

L{r) 



£i{bop{r) + bi[pn{r) - pp{r)]}, 
1 



1 



(III.62) 
(III.63) 

(III.64) 



where A is the Lorentz-Lorenz-Ericson-Ericson correction 
parameter. The kinematical factors €i and £2 are defined 
as £1 = 1 + mT^/rriN and £2 = 1 + mT^/2mN with the nu- 
cleon mass ttiat. The complex parameters Bq and Co are 
the s-wave and p-wave absorption parameters, respec- 
tively. 

From the pionic-atom x-ray data, the p-wave parame- 
ters have been fairly precisely determined, while the x- 
ray data are less sensitive to the s-wave parameters. This 
situation can be understood from Fig. [TOl which shows 
that the Is strong-interaction shift is dominated by the 
s-wave potential while the (and higher) level shifts are 
dominated by the p-wave potential. It is therefore impor- 
tant to obtain experimental information on the Is level 
in order to precisely determine the s-wave parameters. 



2. Pionic hydrogen - the nN scattering lengths at threshold 

In the low-density limit, the potential is described by 
just two parameters, 60 and hi in Eq. pil.62p . which 
are nothing but isoscalar and isovector iiN scattering 
lengths, respectively^. These have been very precisely 



2.5 
2.0 
f 1.5 



m 



1.0 





-0.5 



ii5Sn 



3d(x100) 



Is 



2p(x10) 



FIG. 10 The binding energies with finite-size Coulomb po- 
tential only, Bcoui, and Coulomb plus optical potential, Bfuii, 
are calculated. The energy shifts Bcoui — i^fuii are shown as 
the solid bars for pionic ls,2p, and Sc? states for ^^^Sn. The 
shifts due to the real local terms in the potential are shown by 
dashed bars. Dotted bars are the results wi th all real terms 
(loca l plus nonlocal) in the optical potential (|Umemoto et a/.L 
I2QQQI ). 



determined by the pionic hydrogen x-ray spectr o scop; 
at the Paul Scherrer Institut (PSI) (jGotta et all 1200^ 
ISchroder et al. I. l200lh . 

The ttN scattering lengths can be determined from the 
observed strong-interaction shift ei^ and width Vis of pi- 
onic hydrogen: 



Bis TB 

(X 60-^1, (III.65) 

^ = 8^(1 + 1) k-,_o.(l + M]' 

oc (III.66) 

where = 3.24 keV and = 216 fm respectively are 
the Is binding energy and the 'Bohr radius' of pionic hy- 
drogen, qo = 0.1421 fm~^ the center-of-mass momentum 
of the TT^ in the charge-exchange reaction 7r~p n^n 
and P = 1.546 ± 0.009 the branching ratio of charge 
exchange and radiative capture (Panofsky ratio). The 
quantities Se^r represent the corrections to be applied to 
the experimentally determined scattering length in order 
to obtain pure strong-interaction quantities. 

The PSI experiment R-98.01 (Fig. [11]) used a super- 
conducting cyclotron trap to produce a high stop density 
for pions in a hydrogen gas target, a Johann-type bent- 
crystal spectrometer, and a two-dimensional CCD focal- 
plane detector to determine the pionic- hydrogen x-rays. 
The Is strong-interaction shift was determined from the 
3p - Is energy to be e^f = 7120 ± 11 meV. The Is 



^ lEricson and Ericsonl (|l966l ) defines bo = {ax 12 + 2a32)/3,6i = 
(<^3/2 ~ 0Li/2)/3, while isoscalar a+ and isovector a~ scattering 



lengths are a"^ = (<^i/2 +2 a3^)/3, a~ = ( a-\ — Q-3/2)/3 in pionic 
hydrogen literatures (e.g., (|Gotta et a/ll2008h V Here, ai/2 and 
0^3/2 respectively are the isospin 1/2 and 3/2 scattering lengths. 
Note therefore 60 = a+ and bi = —a~. 
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FIG. 11 A schematic d rawing of the PSI pionic hydrogen 
setup (|Gotta et all l2008l ) . 



hadronic broadening was deduced from the widths of 
Ap — Is^ 3p — Is and 2p — Is transitions, so as to correct 
for the kinematical broadening caused by the preceding 
transitions width, to be Tjf = 823 ± 18 meV. With the 
correcti ons and Sr obtained within chi ral perturbation 
theory (|Gasser et all l2QQ2l : IZemp ^. "2003), the scattering 
lengths have been obtained to be (Got ta et aL , 2008) 



bo = 0.0069 ±0.0031 m-^ , 
bi = -0.0864 ± 0.0012 



(III.67) 
(III.68) 



These values are near the leading-order result derived 
from cu rrent algebra (cahed Tomozawa- W einberg (TW) 
values) (jTomozawai Il966l : IWeinberd . Il966h of 



TW 



0, 



1 TTItt 



-0.079 m: 



(III.69) 
(III.70) 



with ei = 1 + m-j^/rriN = 1.149, m-j^ = 139.57 MeV and 
/tt = 92.4 MeV, revealing an important feature of the 
underlying chiral symmetry. The ~ 10% gap between the 
experimental value of bi and bj^ ^ets closed by pion -loop 



corrections of order ( Bernard et all Il993l , fl995l ). 



3. The missing repulsion problem 

It has long been known that the available pionic-atom 
x-ray data cannot be fitted with bo and bi values, but 
some enhancement of the 5- wave repulsive strength is re- 
quired. This is known as the 'missing repulsion' problem. 

The situation is illustrated in Fig. [121 which the 
5-wave parameters bo and Re^o obtained by various au- 
thors are plotted. We note, firstly, that there is a strong 
correlation between bo and Re^o, approximated by bo + 
0.215 Re^n = —0.028 (known a s the Seki-Masutani re- 
lation (jSeki and Masutanil . [19831 )) ^. Secondly, the Seki- 



^ This linear relation arises since the mean density probed by 
the 7r~ happens to be about 0.5 — 0.6po regardless of element, 



so that bpp -\- KeBpp'^ can be H nearized as 
(|Yamazaki and Hirenzakil . 120031 ) . 
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FIG. 12 The s-wave pa rameters bp vs R eBp ob tained by fit- 
ting the x-ray data; T: (|Tauscher',"l97l'), Bl-3: ("Battx. et all 
1979, 1983, 1997j), SMl-2: (Seki and Masutani, 198j), Kl-2: 
([Koniin et a/.L 11990^). E: pricson and Tauscherl . Il982l i. The 
line is bp + 0.215 ReBp = -0.028. 



Masutani relation implies that the Re^o value must be 
around —0.15 if we are to take the 'free' bo value, which 
is much larger in magnitude than expected from the pion 
deuteron scattering length ( Chatellard et all Il997l ). If 
on the other hand Re^o is forced to be small 0), the 
fitted bo is about —0.03, quite far away from the 'free' 
value. 

Part of the repulsion is known t o be provided by the 
'double scattering' correct ion to bo ([Ericson and EricsonI , 
Il966l : iLoiseau e/oll . l200lh . 



= 



[bl^2bl] 



(III.71) 



where the inverse correlation length is expressed by the 
Fermi momentum kp as 



27r 



3_ 
2^ 



Sir' 



pir) 



1/3 



, (III.72) 



but this correction is still insufficient to account for 
the missing repulsion. The experiment S236 at GSI, 
discussed in the next section, established that the in- 
medium value of bi is enhanced relative to 6i, which, 
through Eq. pil.7ip . provides the extra s-wave repulsion. 



C. Deeply-bound pionic atom spectroscopy 

In 1996, the experiment S160 at GSI reported the first 
observation of th e deeply bound pionic states in ^o^p^ 
( Yamazaki et all Il996h . whose existence and forma- 



tion had__been_Qredj^ 

19851: iHirenzaki et all Il99ll : iToki et all Il991al . Il989 : 
Toki and Yamazaki. " ll988[ ). This discovery opened an en- 
tirely new way to study the hadron properties in the nu- 
clear medium. As illustrated in Fig. [131 the Is wavefunc- 
tion of 7r~ bound to a heavy nucleus overlaps appreciably 
with the nuclear density distribution, and hence the in- 
medium modification of the pion properties may have 
detectable effects on the binding energy and/or width. 
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FIG. 13 a) The pion optical potential for ^°^Pb. The finite- 
size Coulomb potential is expressed by the dotted curve and 
the one with the optical potential by the solid curve. The 
imaginary part is depicted by the dashed curve, b) The pionic 
wavefunctions of the Is, 2s and 2p states in coordinate space. 
The dashed curves and the solid curves are obtained with the 
finite-size Coulomb potential and with the optical potential. 
The half-density radius i?o of ^*^^Pb is indicated by the broken 
line (|Toki ei a/.L Il991bl ) . 



Although this was not clearly recognized when the 
deeply-bou nd pionic-atom spectrosc opy was initially 
conceived ( Toki and Yamazakil , Il988l ), it was later es- 
tablished that the in-medium value of the isovector 
-kN scattering length 6i, derived from the \s bind- 
ing energies, is connected to the (temporal part of 
the) in-medium pion decay constan t ft via the in- 
medium T omozawa- Weinberg relation (Flbmozawal 1 19661 : 
IWeinberd . fT966 ). (Kolomeitsev et a/ll2QQ3bh . 



^i(p) 



47r 



(111.73) 



which is in turn connected to the in-medium chiral con- 
densate via the in- medium Gell-Mann— Oakes— Renner 
(GOR) relation Eqs. (|lL4T]) and (|lL42]) or via the 
Glashow- Weinberg (GW) relation Eq. (|IL431) . 

In the case of deeply-bound pionic atoms, the i\~ quan- 
tum number is well defined. Moreover, as will be shown, 
the use of recoilless {d^ He) kinematics ensures that the 
nucleus is in the ground state. The pion wavefunction 
and its overlap with the nuclear density are therefore pre- 
cisely calculable by solving the Klein-Gordon equation. 
This makes it possible to compare the experimental re- 
sults with theoretical predictions, and to quantitatively 
deduce the effects of partial restoration of chiral symme- 
try. 




Deeply bound region 



> r^bs > r^ad 



FIG. 14 A schematic figure of pionic atom states with x-ray 
transitions down to the last orbital and deeply bound inner 
orbits with large widths which cannot be populated following 
the x-ray cascade. They have large widths due to nuclear 
absorption but a r e still discrete states with Fn < En — En-i- 
(jYamazaki et a/.Ll2008l ). 



A wealth of data on pionic atoms has b een collected 
by m eans of pionic x-ray spectroscopy ( Battv et all 
1997r). In the x-r ay spectroscopy experiments (e.^ 



de Laat et al pions, injected into a target, slow 

down and form pionic atoms, and x-rays emitted in 
the cascade are measured, as schematically depicted in 
Fig. [HI While high-lying states do not show detectable 
strong-interaction effects, the energy levels of low- lying 
states get shifted from those calculated by using the elec- 
tromagnetic interaction. The level widths also become 
larger due to pion absorption on the nucleus. From the 
measured strong-interaction shifts (e) and widths (F), the 
pion-nucleus strong-interaction potential parameters can 
be deduced. 

The measured pion-nucleus strong-interaction shifts 
are shown in Fig. [15] as function of atomic number Z. 
This plot shows that the information on the Is level shift 
is available only up to the atomic number of Z = 14. This 
is because beyond the last orbital (Fig. [T4j) x-rays can- 
not be observed due to the increase of the pion-nucleus 
absorption width. As shown in Fig. [151 3(i is the last or- 
bital for medium-to-heavy nuclei, meaning that it is not 
possible to obtain the information on Is and 2p levels by 
means of x-ray spectroscopy. The Is and 2p levels which 
cannot be populated by x-ray cascade are often referred 
to as the deeply-bound pionic states. 



1. Structure of deeply-bound pionic atoms 

This, however, does not mean that those states do not 
ex ist, as was rec ognized bv Friedman and Sofj (|l985l ) and 
bv lToki and Ya mazaki (19881). In Fig.[T6l a tvpical pion- 
nucleus optical potential set (see section IIII.Bj) was used 
to calculate the binding energies and widths of pionic Zr 
(left) and Pb (right) atoms. The Is level width Fi^ is 
found to be much smaller than the 2p — Is level interval 
{AE2p-is < Tis) even for heavy pionic atoms such as 
pionic ^^^Pb. Namely, the deeply-bound pionic atoms are 
metastable, despite the strong pion-nucleus absorption. 

This somewhat counter-intuitive result arises due to 
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FIG. 15 Strong- interaction shift (top) and width (bottom) 
of pionic atoms. The continuous l ines join points ca lculated 
with the best-fit optical potential (|Battv et all Il997l ) . 



the repulsive nature of the pion-nucleus s-wave interac- 
tion, which causes the Is binding energies to decrease 
(Fig. [T6|) and also pushes the wavefunction outwards (see 
the difference between the solid and dashed curves in 
Fig. [T3lf b)). This then reduces the pion-nucleus overlap, 
making the level width narrower. 



2. Formation of deeply-bound pionic atoms 

The original dee ply-bound pionic atom fo rmation 
scheme proposed by iToki and Yamazakil ( 19881 ) was to 
use the ^Q^Pbfn p) re action. This was tested at TRIUMF 
(llwasaki et a/.LIl99lh . but no bound-state peak was ob- 
served. It was soon recognized that (n, d) or (d,^ He) re- 
actio ns (see Fig. [T7|) are more suitable (|Hirenzaki et all 
Il99l[ ). These are recoilless as well as substitutional re- 
actions, in which a neutron in the s shell is picked up 
and the produced 7t~ is left in the Is orbit with a small 
momentum transfer (Fig. [T8|). Due to the substitutional 
nature of the reaction, the angular momentum transfer is 
AL = 0. The small momentum transfer makes it possi- 
ble to satisfy the angular momentum matching in surface 
reactions, thereby minimizing the nuclear distortion. 

Using the (d,^ He) reaction, the experiment S160 at 
GSI succeeded to observe the 2p and Is states of pionic 
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FIG. 16 Energy levels of pionic atoms of ^°Zr and ^°^Pb. The 
results with the finite-size Coulomb potential (i.e., taking ac- 
count of the nuclear charge distribution) are shown by dashed 
bars, while those with the optical potent ial by solid bars with 
hatched area indicating the level widths (|Toki and Yamazakil . 
119881 ). 
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FIG. 17 Diagrams for proton-pick-up pion-transfer (a) (n, d) 
and (b) ((i,^He) re actions to form pionic bound states on a 
neutron-hole state (|Hirenzaki et all Il99ll ) . 



2°^Pb (jYamazaki et a/.LIl996h . As shown in Fig. [IS the 

experimental result agrees remarkably well with the the- 
oretical prediction. This firmly established the method- 
ology of deeply-bound pionic atom spectroscopy. 



3. GSI S236 - Sn(d,2He) 

In principle, the 5- wave parameters may be obtained 
by analyzing the Is-level energies and widths of light pi- 
onic atom x-rays {Z < 14, Fig. [T5|). However, since the 
available data are on TV = Z nuclei, the sensitivity of the 
dataset to the isovector parameter bi is quite limited. 

The S236 experiment at GSI, a successor to S160, 
measured the Is binding energies and widths of pionic 
ii5,ii9,i23g^ using the ii6'i20,i24Sn(^^3He) reactions. Sn 
isotopes were chosen since the Is states are expected to 
be produced as the most dominant quasi-substitutional 
states, (15)7^- (3s)~^, because of the presence of the 35 
orbital near the Fermi surface (Fig. [2Q|). Another merit 
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FIG. 18 Momentum transfers in the (n,d) and (d,'^He) reac- 
tions on Pb for Q — —130 and —140 MeV as a functi ons of 
the incident energy per nucleon (|Hirenzaki et alV Il99lh . 
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FIG. 19 (a) The ^°^Pb(d,^ He) spectrum. The abscissa is 
the reaction Q value, and the tt" emission threshold {Q — 
-140.14Mey) is shown by a dotted line. The He)7r° 
reaction peak used as a calibration is also shown in the Q value 
scale of the ^° ^Pb((i,^ He) reac t ion ki nematics, (b) Theoretical 
prediction bv lHirenzaki et a/ 1 ([T99ll) . A F^y HM resolution of 
0.5 MeV is assumed (|Yamazaki et a/.LIl996l ). 



is to make use of isotopes over a wide range oi{N — Z)/A 
to sensitively deduce the isovector parameter hi. 

The experiment S236 used a deuteron beam from the 
heavy- ion synchrotron SIS at GSI, Darmstadt, combined 
with the fragment separator (FRS) as a high-resolution 
forward spectrometer (Fig.[2T|). They chose the exact re- 
coilless condition to suppress minor states other than the 
enhanced Is states with quasi-substitutional Is states, 
with a deuteron beam of a small momentum spread and 
an accurately measured energy of 503.388 ± 0.100 MeV. 
The Q-value resolution was 394 ± 33 keV (FWHM), and 
the absolute Q-value scale was calibrated to an accuracy 
of ±7 keV using the p((i,^He)7r^ reaction, where a thin 
Mylar layer put on the surface of each Sn target was used 
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FIG. 20 The Sn isotope dependence of the total (d,^ He) spec- 
trum for the pionic atom formation at Td — 500 MeV with 
300 keV experiment al resolution. The targ et nucleus is indi- 
cated in the figure. (|Umemoto et all . 12000 
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FIG. 21 The Fragment Separator of GSI used for the (d,^He) 
spectroscopy. 



as the proton source. 

The observed spectra, j dE j dVt^ on Mylar-covered 
"'^"'^^Sn, -"^^^Sn, -'^^^Sn targets as a function of the ^He ki- 
netic energy, are shown in Fig. [22l In each spectrum 
of Fig. [22] a distinct peak at around 365 MeV was ob- 
served, which was assigned to a dominant configuration 
of (ls)7r(3s)~^. The skewed peaks at around 371 MeV 
arise from p((i,^He)7r^. The overall spectrum shapes for 
the three Sn targets were found to be in good agree- 
ment with the predicted ones (Fig.[20|). The spectra were 
decomposed according to the theoretical prescription of 
lUmemoto et al\ ( 2QQQ| ). from which the Is binding ener- 
gies [B\s) and widths (Fis) were determined as: 



3.906 ± 0.021 (stat) ± 0.012 (syst) MeV, 
0.441 ± 0.068 (stat) ± 0.054 (syst) MeV, 
3.820 ± 0.013 (stat) ± 0.012 (syst) MeV, 
0.326 ± 0.047 (stat) ± 0.065 (syst) MeV 
3.744 ± 0.013 (stat) ± 0.012 (syst) MeV, 
0.341 ± 0.036 (stat) ± 0.063 (syst) MeV. 
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FIG. 22 Double differential cross sections versus the "^He ki- 
netic energy of the ^^'^'^^°'^^^Sn((i,^He) reactions measured at 
the incident deuteron energy of 503.388 MeV. The scales of 
the 7r ~ binding energies in 123,119,115 ^^^^ indicated. 
From ISuzuki et all (|2004l l . 



In-medium isovector scattering length bi 

The s-wave potential parameters {60, 61, Re^o, Im^o} 
were deduced by simultaneously fitting Bis and Tig of the 
three Sn isotopes together with those of symmetric light 
nuclei (^^O, ^^Ne and ^^Si), with the p-wave parameters 
fi xed to the known v alues from pionic x-ray data listed 
in lBattv et all (|l997h . 

The obtained values are: 



^0 

bi 

ReBo 
ImBo 



-0.0233 ±0.0038 m-^ , 
-0.1149 ±0.0074 , 
-0.019 ±0.017 m-^ , 
0.0472 ± 0.0013 m-^ . 



The effective density probed by the 7t~ was p ^ 
0.6po (obtained from the overlap of the nuclear den- 
sity and the pion wavefunction , such as shown in 
Fig. [13] (|Yamazaki and Hirenzakil . l2003h ). Figure [25] 
shows the likelihood contours in the plane of {61, Im^o}- 
As shown, the best-fit bi value deviates significantly from 
the 'free' value obtained from pionic hydrogen. 

Also shown therein is how the best-fit values move if 
the shape of the neutron distribution was changed. In 
fact, the most serious problem in the analysis was the 
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FIG. 23 Difference Arnp between the rms radii of the neu- 
tron and proton distributions as deduced from the antipro- 
tonic atom x-ray data, as a function of ^ = (N — Z)/A. The 
proto n distributions were obtained from electron scattering 
data ( De Vries et a/.LIl987l ) (Sn nuclei) or from muonic atom 
data ( Fricke et ami995l l (other nuclei). The full line repre- 



sents the linear relationship between 6 and Avnv as obtained 
from a fit to the experimental data. From iTrzcinska et all 
(|200lh . 
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FIG. 24 Difference Arnp between the rms radii of the neutron 
and proton distributions for various ti n isotopes deduced from 
proton elastic scattering at 295 MeV (|Terashima et all l2008l ) 
(s quares) comp ared with the values adopted in the analysis 
of lSuzuki et al.\ {[2004) (filled circles), which used the relation 
Arnv = (1.01 =b 15)( A^-Z)M+ (-0.04 zb 0.03) obtained by 
ITrzcinska et al\ (|2Q0lh . 



relatively poor knowledge concerning the neutron distri- 
bution Pn{r) in Sn isotop es, whereas the prot on distribu- 
tion pp{r) is well known ([Fricke et all Il995l ). In the fit, 
the neutron diffuseness (a^) and half-density radius (c^) 
parameters in the two-parameter Fermi distribution were 
chosen so as to satisfy the difference between the neutron 
and proton rms radii, Arnp = (1-01 ± 0.15){N - Z)/A + 
(— 0.04±0.03) fm (Fig.[23|). based on experimental data of 
antiprotonic atoms of Sn isotopes (Trzcihska et al.^ 2QQlh 
as well as of many other nuclei. The two extreme assump- 
tions, the "skin" type (cp < Cn^dp = an) and the "halo" 
type {cp = Ca^ap < an)-, where Cp (cn) is the proton 
(neutron) half-density radius and ap (an) is the proton 
(neutron) diffuseness of the two-parameter Fermi distri- 
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FIG. 25 Likelihood contours in the {bi,lmBo} plane from 
the simultaneous fitting of {Bis^Tis} of t he Is pionic states 
in the three Sn iso topes and three light symmetric nuclei 
([Suzuki all I2QQ4I ). 



bution, are shown in Fig. [25] in open diamonds and open 
squares, respectively. The best-fit values were obtained 
by using the parameters halfway between these two as- 
sumptions. The two crosses in the figure indicate the 
dependence of bi on the uncertainty of ±0.04 fm in A^p. 
The quoted errors on the best-fit values do not include 
these uncertainties in the neutron distribution^. 



In-medium quark condensate 

From the observed enhancement of the hi parame- 
ter relative to the free value, bi/bi{p) = 0.78 ± 0.05, 
ISuzuki et al. 1 (I2OO4D deduced that the chiral order pa- 
rameter is subject to the in- medium reduction of 
{fl{p)/ f^)"^ 0.64 at the n ormal nuclear density p = po , 
based on th e sugg estion ( Kienle and Yamazakil , I2OOII : 
IWeisd . [2000[ l2001bh that the missing repulsion may be 
explained in terms of a possible in- medium change of t he 
pion decay constant. A global fit ( Friedmai] , l2002al lbl) 



to th e pionic x-ray as w ell as the deeply-bound ^^^Pb 
data ( Geissel et all |2002[ ) supported this view. Recent 



direct calculation of 61/61 (p) in the unitarized chiral ap- 
proach (Fig. [261) is also consistent with the above result 
(jPoring and^ Oset. 2008). 

Here, we adopt a new model independent rela- 
tion obtained by exploiting operator relations in QCD, 
Eg. (1X1431) . IJido et al\ (|2008[ ) showed that Zl^^{p) 

th 7 = 0.184. Using 61/61 (p) = 0.78 at 

p = 0.6po7 the ratio of the quark condensates is found to 
be 



1-0.37- 



Po 



(III. 74) 



^ Recently, iTerashima et al\ (120081 ) have determined the neutron 
rms radii of Sn isotopes using proton elastic scattering at 295 
Me V. Their resu l ts are in good agreement with the radii assumed 
in (ISuzuki et all\2004 ). See Fig. [24] 
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FIG. 26 In-medium isovector bi{p) compared to the vacuum 
isove ctor term bi. The gray band, the experimental result 
from [Suzuki et q H ("2004'), is co mpared with chiral calcula- 
tions by IPoring a nd Oseti ([20081) (^dark ban d ) . [Meissner et al\ 
(I2OO2I ) (dashed hne) and I Friedman et al\ (|2Q0 4|) (^dash-dot 
line) . The p oint is a phenoni e nolo^ ical fit bv iNieves et al\ 
(I1993I ). From lDoring and Osetl (l2008l ) 



A similar conclusion, (qq) ^ / (qq) — 1 — 0.39/)//)o, was ob- 
tained from the differential cross sections for tt^ -nucleus 



elastic scattering data at 21.5 MeV ([Friedman et 
[2005. .2004 ). 



IV. SCALAR MESON: a IN NUCLEI 
A. Theoretical background 

The <j-meson in the vacuum is a very broad resonance 
in the scalar-isoscalar channel as discussed in Sec lII.B.2l 
It is an excitation having the same quantum number as 
the vacuum, and may be interpreted as the Higgs boson 
in QCD. The fate of such scalar-isoscalar excitation in 
hot and/or dense medium i s strongly correlated wi t h chi- 
ral symmetry restoration ( Hatsuda and Kunihird [l985[ , 
1987a). The basic idea is rather simple and general: 
If a sizable reduction of the chiral condensate takes 
place in the medium, the ground state becomes "soft" 
against the amplitude fluctuation of the order parame- 
ter due to the reduction of the stiffness. In fact, the 
frequency of the a excitation (the amplitude fluctua- 
tion) is red-shifted toward the frequency of the tt ex- 
citation (the phase fluctuation). Similar softening phe- 
nomena are well known in solid state physics, e.g., the 
soft pho non modes in anti ferro-el astic crystals such as 
SrTiOs ([Gebhardt and Krev. "l980) and in ferro-electric 
crystals such as SbSI (Kittel, 2004). 

In the case of QCD, an interesting signal associ- 
ated with the chiral softening is the spectral enhance- 
ment in the scalar-isoscalar channel near the 27r thresh- 
old at finite temperature and density as mentioned in 
Sec lII.C.Tl Consider the retarded propagator of a at 
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rest in the medium, D(j{uj). The spectral function is 
defined as = -^ImD^ = ^ImD'^ /[(ReD-^)'^ + 
{lmD~^)'^]. Because of the softening and strong airir cou- 
phng, IleD~^{cj 2mT^) becomes smah or even vanishes 
at a certain temperature or baryon density. In that situa- 
tion, the spectral function is dominated by the imaginary 
part of the invers e propagator wit h the phase space fac- 
tor (IChiku and Hatsuda . 1998a .b: iHatsuda eta/] . Il999l : 
IVolkov et a/.l . ll998[ ). 
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This implies a large enhancement of the spectral function 
near 27r threshold. Such an enhancement may be seen in, 
e.g., the dipion production and diphoton production from 
the hot/dense medium. Also, a a-mesic nuclei (or the 
bound dipion in nuclei) could be formed by {d^ t), {d^ ^He) 
and (7,p) reactions, if there is large enoug h softening 
( Hirenzaki et l2QQ2l : In asahiro et al. I l2QQ5h . 

In reality, the pion has a width inside the medium, 
so that the spectral function in the a-channel does not 
have a simple form such as Eq. pV.75p . Nevertheless, 
more sophisticated approaches indicate similar enhance- 
ment of the in-medium tttt scattering amplitude in the 
scalar-isoscalar channel near the 2tt threshold (see the 
references cited in Sec JII.C.7|) . Shown in Fig. [271 is one 
of such examples obtained by using the chiral unitary 
model. Here we note that the expansion parameter of the 
non-linear chiral models is a ratio, [typical momentum 
of the pion or nucleon fields]/ [chiral symmetry break- 
ing scale 4:7Tf-j^]. Therefore, non-linear chiral approaches 
loose their predictive power at high temperature and/or 
high baryon density. 



B. Experiments 

1. CHAOS 

The CHAOS (Canadian high acceptance orbit spec- 
trometer) collaboration at TRIUMF was the first to re- 
port such a low- mass enhancement in the (7r7r)/=j=o 
channel. They used 7r~^A. ^ tt^tt~ X reactions on hy- 
drogen (|Kerm ani et al.. Il998l) and on nuclear targets 
(iBonutti et a l. 199 6. Il997l Il999l . |2QQQ| : ICamerini etal\ , 
l2QQ4i l2QQll : TGrion et ali l2QQ5h at pion kinetic energy 
= 243 - 305 MeV, and observed the spectral "soft- 
ening" in the tt+tt" channel but not in the 7r+7r+ (i.e., 
1 = 2) channel. The nuclear data were taken at = 
283 MeV on ^H, ^^c, ^OCa and ^o^Pb targets, as weh as 
at = 243, 264, 284, 305 MeV on ^^Sc. 

Figure [28] shows a typical reconstructed tt+tt" event. 
The spectrometer is based on a cylindrical dipole mag- 
net producing vertical magnetic fields up to 1.6 T (0.5 T 
for the 27r experiments) ([Smith et a/.LIl995h . The target 
is located in the center of the magnet. Charged particle 
tracks produced by pion interactions are identified us- 
ing four concentric cylindrical wire chambers (WCl, 2, 
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FIG. 27 Imaginary part of the tt — tt amplitude Tqq{E) in the 
scalar- isoscalar channel calculated in the chiral unitary model 
(jCabrera et ali E0O8) . 
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FIG. 28 A typical reconstructed event in CHAOS for the 
Tif 7r~^7T~p reaction on ^^C. WCs are the four wire cham- 
bers, CFTs are the CHAOS First-level Trigger counters, 
which are also used for particle identification. Of the 20 CFTs, 
two a re removed to free the pion beam path. (|Bonutti et all 
I2QQQI ). 



3, 4) surrounding the target. Particles are identified by 
cylindrical layers of scintillation counters and lead-glass 
Cerenkov counters, which also provide a first level trig- 
ger (CFTs). The detector subtends approximately 10% 
of 47r. The momentum resolution delivered by the detec- 
tor system is 1%. The pion-detection threshold energy is 
11 MeV. 

Fig. [29] shows the tt+tt" (left) and 7r+7r+ (right) 
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FIG. 29 Invariant mass distributions (diamonds) for the 
TT^ 7i~^7i~ and TT^ TT^TT^ reactions on ^H, ""^^C, ^°Ca and 
^°^Pb. The shaded regions represent the results of phase- 
space simulation s for the pion-p roduction reaction ttA 
7T7tN[A - 1] (jBonutti et a/J . 120001 ). 



invariant- mass spectra taken on nuclear targets. The dis- 
tributions span the range from 2m-j^ up to the 420 MeV, 
the maximum allowed by the reaction. While the 7r+7r+ 
spectra can be fairly well represented by the phase-space 
simulations (shaded region), the tt+tt" spectra show a 
peak in the low-mass region, increasing with mass num- 
ber A. 

Instead of comparing the raw spectra with theoretical 
predictions or with results of other experiments, a com- 
posite observable was used: 
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FIG. 30 (left): C:^^, the bin- by-bin ratio (see text) of tttt 
invariant mass distributions for the two reactions tv^A 
-K^-K^ A! and vr^p Tr^vr^n, as a function of the Mtttt energy. 
The nuclei (A) examined are (which plays the role of a 
proton, p), ^^C (open diamonds), ^°Ca (full triangles), and 
^°^Pb (open squares), (right): Sa me as the left panel b ut for 
the TT^ tt^tt" reaction channel (|Bonutti ei a/.L[l999l ). 



2. Crystal Ball 

The Crystal Ball (CB) collaboration at the AGS in- 
vestigated the 7r~ A ^ tt^tt^A' reaction on CH2, CD2, C , 
Al an d Cu targets at p^^- = 408 MeV/c (|Starostin et all 
|2QQQ[ ). The Crystal Ball comprises 672 optically isolated 
Nal(Tl) crystals that cover 93% of 47r, and has a tt^tt^ 
acceptance of 11 — 17%. The invariant mass resolution 

is about 1.2% at m^o^o = 2m^o and reaches a plateau 

0^0 



of 2.2% at 



0.3 GeV/c^. The measured tt^tt' 



invariant- mass spectra (left panel of Fig.[3T]) show a grad- 
ual shift of intensity towards lower m-j^-j^ for heavier tar- 
gets, but a sharp strong peak near 2m^ as reported by 
the CHAO S collaboration cannot be seen. 

However , Gamer ini et al\ (|200lh pointed out that if the 
composite ratios are used to compare the CHAOS 
and CB results, so as to (mostly) remove uncertainties 
arising from acceptance corrections^, the two results are 
not statistically inconsistent, at least in the case of -^^C 
(the only nucleus common to the two experiments). See 
Fig. [31] right panel. 



3. TAPS 



where cr(M^^) ((j(M^)) is the triple differential cross 
section d^a/dM^^drt^d^}^^ for nuclei (nucleon), M^^ rep- 
resents the TTTT invariant mass, denotes the pion- 
detection solid angle, and {(^t) total cross 

section in nuclei (nucleon). The ratios are presented 
in Fig. [30l The ratios C^+^- (right panel) show that 
the low-mass tt+tt" pairs are more abundant in heavier 
targets, while no such trend can be seen in C^+^+ (left 
panel). 



The TAPS (Two Arms Photon Spectrometer) collabo- 
ration us ed the tagged photon facility at the MAMI ac- 
celerator ( Anthonv a/.lJl99ll ) to me asure A (7, tt^tt^) as 
well as ^4(7, TT^TT^/") cros s sections (|Bloch et al 1 12OO7I : 
iMesschendorp et a/.L [20021 ). 

The experiment covered the photon energy range from 
200800 MeV with an energy resolution of 2 MeV per 



^ Note that the CHAOS acceptance is about 10% of An while that 
of Crystal Ball is 93% of 47r. The small acceptance of CHAOS 
may be the origin of the sharp peaks close to the threshold. 
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FIG. 31 (Left) Experimental results for the 27r° invariant 
mass distributions obtained for the H, D, C, Al, and Cu 
targets corrected for Crystal Ball acceptance. The verti- 
cal scale is in arbitrary units. The solid lines show the re- 
sults of calculations made for the TAPS group by Rapp, 
and the dashed line is the prediction by Vicente Vacas. 
(jStarostin et all l2000l l . (Right) The composite ratios C:^^ 
as a function of the tttt invariant mass for the ^^C target. Full 
di amonds, the C^o^o di stribution deduced from the CB data 
of lStarostin et aU (|2QQQ | ): open diamonds, th e CHAOS 
distribution taken from Bonutti et al.\ (| 19991 1 ([Camerini et all 
[2OOI1). 



ta^^er channel. The targe ts used were ^H , ^^C, "^^Pb 
(iMesschendorp et all\20m as well as ^^Ca ( IBloch et all 
I2QQ7D . The reaction products from the target were de- 
tected with the electromagnetic calorimeter TAPS, com- 
prising 510 hexagonally shaped BaF2 crystals of 25 cm 
length with an inner diameter of 5.9 cm. They were 
arranged in six blocks of 64 modules and a larger for- 
ward wall of 138 modules (see Fig. [32|) . The blocks 
were arranged in one plane around the target at a dis- 
tance of 55 cm from the target center and at polar angles 
of ±54°, ±103°, and ±153°, while the forward wall was 
placed 60 cm away from the target center at 0° and the 
photon beam passed through a hole in the center of the 
forward wall. Each detector module was equipped with 
an individual plastic veto detector, read out by a sepa- 
rate photomultiplier. The setup covered 37% of the full 
solid angle. The two-7r^ invariant- mass resolution was 
between 2.0% and 2.5% in the incident-photon energy 
range of interest. 

Figure [33] shows the tt^tt^ (left) and tt^tt+Z- (right) 
spectra measured at the incident photon energy range of 
400 — 460 MeV. This energy range was chosen so that 
its centroid corresponds to the same center-of-mass en- 
ergy as was used in the pion-induced experiments, en- 
abling a direct comparison. Since this range is below the 
/^-production threshold of 550 MeV, the event identifica- 
tion is clean. Figure [33] indicates that the strength in 
the distribution of M^o^o (but not M^o^+/-) is shifted 



towards smaller invariant masses with increasing A. 

A more r ecent analys is of ^ ^Ca(7, tttt), with higher 
statistics by iBloch et all (|2007h . shown in Fig. [34l re- 
vealed that the invariant-mass spectra show a similar 




FIG. 32 Setup of the TAPS detector at the Mainz MAMI 
accelerator. The beam entered the target chamber f rom t he 
lower right edge. (|Bloch et all l2007l : iKermani efall 1 19981 ) 
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FIG. 33 (Left) Mass- number- normalized differential cross 
sections of the reaction A(7,7r°7r°) with A = ^H, ^^C, ""^^Pb. 
(Right) The same for the reaction A(7,7r°7r^). Both panels 
are for incident photons in the energy range of 400 — 460 MeV 
(soUd circles). Error bars denote statistical uncertainties. 
The dotted curves indicated phase-space distributions deter- 
mined by the Monte Carlo mod el, and the solid curves are 
predictions by (|Roca et all\200'i ). From [Messchendorp et all 
(,2002 ). 



soften ing effect as already found in [Messchendorp et all 
( 2002[ ) for carbon and lead nuclei ^, and that the strength 
of the effect is comparable to carbon, but they also found 
that a sizable part of the in-medium effects can be ex- 
plained by final-state interaction effects, namely, pion 
rescattering, as discussed below. 



^ Note that the cross-section ratios f7(7r*-'7r^)/cr(7r°7r*-') are about 
2 in Fig. [33] while they are about 5 in Fig. [34] This is likely 
due to the larger sys t emat ic errors in the tt ^tt^ cross sections in 
IMesschendorp et al\ (|2002l ) (|Kruschd . [Soosi ) . 
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FIG. 34 Pion-pion invaria nt- mass distribut ions compared to 
results of the BUU model (iBuss et a/.Ll2006l ). The bars at the 
bottom represent the systematic uncertainty of the data, the 
dashed lines r epres ent the error band for the BUU calculation 
(|Bloch et ad 120071 ). 




FIG. 35 (color) Two-pion invariant mass distributions for 
7r°7r° photoproduction off ^^C and ^°^Pb for = 0.4 - 
0.46 GeV, calculated bv iBuss oil (|2006l ). Results without 
final-state interactions (black dash-dotted lines), and with FSI 
(black solid lines). These are compared with the results ob- 
tained by (Roca et a/., 2002), with in- medium final tvtv inter- 
action (red solid lines) and with the free-space tttt interaction 
(red dashed lines). 



C. Final-state interaction (FSI) effects 



The solid cur ves in Fig. [33] are predictions by 
iRoca et a/] (|2002[ ). Here, the meson- meson interaction 
in the I = J = channel is studied in the framework of 
a chiral-unitary approach at finite baryon density. The 
model dynamically generates the a resonance, reproduc- 
ing the meson-meson phase shifts in vacuum and ac- 
counts for the absorption of the pions in the nucleus. 
In the model, the tt — tt FSI modified by the nuclear 
medium produces a shift of strength of the tttt invariant 
mass distribution induced by the moving of the a pole to 
lower masses and widths as the nuclear density increases. 
The data are described well by the model considering a 
theoretical uncertainty of 20%. 



On the other hand, the curves in Fig. [3ll are the re- 
sults of the semi-cla ssical Boltzmann- Uehling-Uhlenbeck 
(BUU) calculation (|Buss et all [2006'), which reproduce 
both TT^TT^ and tt^tt+Z" data reasonably well. This model 
does not contain the tt — tt final-state interactions, but 
the "softening" of the tt^tt^ spectra is due to charge- 
exchange pion-nucleon scattering (i.e., tt — N FSI) which 
mixes the contributions from the different charge chan- 
nels. Since the total cross-section for tt^tt^ production 
is much larger than the tt^tt^ cross-section, the latter 
receives significant side feeding from the mixed charge 
channel via tt^N ^ tt^N scattering, which increases the 
fraction of re-scattered low-energy pions in this channel. 
In the same way, re-scattering of tt+tt" contributes to 
the TT^TT^ channel. 



V. VECTOR MESONS: p,cj,0 IN NULCEI 
A. Theoretical background 



As is discussed in Sec JII.B.3l the direct signature of chi- 
ral restoration is the degeneracy between the vector spec- 
tral function and axial- vector spectral function p^. 
Since the vector current couples to virtual photons which 
eventually decay into dileptons (/^/~), Py is directly re- 
lated to the physical observable. For example, the emis- 
sion rate of dileptons (number of dileptons emitted per 
space-time volume d'^x and per energy-momentum vol- 
ume d^p) from the hot /dense matter reads 



d^xd'^p SttV e'^/T - 1 



I{mf/p^). (V.76) 



Here = (cc;,p) is the total four momentum of 
and the superscript T (L) implies transverse (lon- 
gitudinal) and I{z) = (1 + 2z){l - ^zf 1^0(1 - 4z) with 
z = mf/p'^ denotes the phase space correction from the 
finite lepton mass, mi. Unlike the case of p^, it is dif- 
ficult to measure p^ by dilepton, since the decay occurs 
through and is highly suppressed at low energies. 

The spectral constraints on p^ it self are obtained from 
the operator product e xpansion ([Hatsuda et al I [1991; 
iHatsuda and Led , Il992l ) similar to the derivation of the 
Weinberg- type sum rules in Sec lII.B.3l 

^(PvM-PpqcbM)=0, (V.77) 

poo 

/ du;2(p,(c.)-p^„,,H) = ^Cl{01), (V.78) 

-'0 
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poo 

/ du;W{pAu^) - Ppcco(^)) =J2^'M). (V.79) 



Here (9^ are the local composite operators with dimen- 
sion n with Lorentz indices in general, and are the cor- 
responding Wilson coefficients. In the SU(2) chiral limit 
^u,d = are all chirally symmetric, while 0\ contain 

both chirally symmetric and non-symmetric operators. 
Also, Ppqcd(^) is the spectral function which reproduces 
the perturbative calculation of the correlation function 
(the l.h.s. of Eq. pi.3ip ) in the deep Euclidean region 
(cj^ —00). Then, Ppqcd(^) chirally symmetric by 
definition. These are the reasons why the Weinberg-type 
sum rules with only chirally asymmetric condensates in 
Sec JII.B.3l are obtained by taking the difference between 
vector and axial- vector correlations. Since chirally sym- 
metric condensates do not have to vanish at the critical 
point of the chiral transition, one cannot immediately 
relate the spectral modification of vector mesons to the 
restoration of chiral symmetry. 

Even if in- medium changes of {0\q) are obtained ex- 
actly from lattice QCD simulations, the above sum rules 
only supply information on the weighted averages of the 
spectral function and not on the exact spectral shape. 
Nevertheless, these sum rules are useful to make a con- 
sistency check in various models of QCD (jKlingl et all 
I1997I : iKwon e^"all , l2QQ8l ). Also, these sum rules may be 
used to extract the information on {0\q) by adopting the 
experimental dilepton spectrum after background sub - 
traction in the l.h.s. of Eqs. (|V.78IV.79P (jHatsudal . [TqQTI ) . 

In general, the spectral function receives peak-shift, 
broadening, new peaks, etc due to the complex in- 
teraction of the vector current with the medium 
( Rapp and WambachlJ2QQQl ). Also, such spectral changes 
may wel l depend o n the spati al momen t um o f the current 
(^Elgtsk v and Ioffel .fT997: Frim an et a/lll999l : [jean et all 
I1994I ). Indeed, the transverse and longitudinal spec- 
tral functions p^'^(cj,p) obey different Weinberg-type 
sum rules for p ^ (|Kapusta and Shurvakl . [T99l . Be- 
cause of these reasons, it is not appropriate to oversim- 
plify the problem to "mass shift vs. width broaden- 
ing" . Nevertheless, there is a theoretical suggestion that 
the width broadening at low temperature and/or baryon 
density is eventually taken over by th e mass shift near 
the critical point of chira l transition (|Yokokawa et all 
l2QQ2[ ) (|Brown et all l2QQ8h . Experimentally, it is im- 
portant to measure the full momentum dependence of 
the spectral function p^'^(c(;,p) instead of the projected 
invariant- mass spectra. 



ity of decaying in medium, while co and (j) tend to decay 
outside. In order to study co/cj) in- medium properties, it 
is important to choose a proper reaction and to select 
slow- moving mesons. 

Obtaining dilepton distributions (usually presented in 
the form of l~^l~ invariant- mass spectra M1+1-) is techni- 
cally very demanding because of the small dilepton-decay 
branching ratios of these mesons (Tab. [Tj), while there 
are many hadronic sources which can produce leptons. 
The detector therefore must have an excellent lepton- 
identification capability, and must also provide means 
to suppress combinatorial background, the background 
caused by an being erroneously paired up with an /~ 
from other origin (e.g., a e+ from tt^ 76+ e~ paired up 
with an e~ from 7 ^ e+e~ occurring in the same event). 

Even with the start-of-the-art dilepton detectors, the 
combinatorial background is severe, especially in high- 
energy heavy- ion collisions. For example, at CERN 
SPS (158 AGeV central cohisions), the NA60 exper- 
iment (/i+/i~) report ed a signal S to b ackground B 
ratio of about 1/11 ([Damianovid . l2QQ7h . and it was 
about 1/22 in th e case of the CERES experiment (e+e~) 
(lAdamova et all I2QQ8D . At RHIC (y^ = 200 GeV 



B. Dileptons, why and how? 

Dileptons (/+/~ pairs, where / = e or /i) are an ex- 
cellent tool to study possible in-medium modifications of 
vector mesons (Table H]), p^co and (j) in nuclear media, 
because of their negligible final-state interactions. Due 
to their short lifetime, mesons have larger probabil- 



Au-hAu minimum bias), the PHENIX experiment (e+e~) 
reported a signal-to-b ackground ratio of about 1/100 
(|Afanasiev et all\200% . 

In order to reliably extract meaningful results despite 
such small signal-to-background ratios, methods such as 
event mixing and like-sign pair subtraction have been de- 
veloped to reliably subtract combinatorics, as discussed 
in section [V.B.3[ 

The combinatorics-subtracted Mi+i- distribution still 
contains a broad continuous background due to Dalitz 
decays. In order to extract the vector- meson contribu- 
tions, the measured distribution is compared with the 
"hadronic cocktail" , which contains all known sources of 
pairs produced in the detector acceptance. 



1. /i^/i -pair detection 

Muons produced in pion and kaon decays (e.g., 7r+ 
p^v and M^^) are much more abundant than 

those from vector-meson decays, and they contribute to 
the combinatorial background. It is therefore essential 
to absorb hadrons as close as possible to the interaction 
point in a thick absorber. Muon momenta are measured 
by magnetic spectrometer(s) placed behind the absorber, 
corrected for the energy loss and multiple scattering in 
the absorber, and the pair-mass distribution is recon- 
str ucted. This is exa ctly how T{hh) was first discovered 
bv lHerb et al\ (|l977D . 

This method works well for heavy-mass region (T as 
well as J/^), where the decay-muons have high mo- 
menta, and the particle multiplicities are low, but the 
measurement of low-mass vector mesons is difficult due 
to larger combinatorial background and larger multiple 
scattering in the absorber (hence lower mass resolution). 
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TABLE I Properties of vector mesons. (|Amsler et all I2QQ8I ) 





Mass 
(MeV/c^) 


r 

(MeV/c^) 


CT 

(fm) 


Main 
decay 


Ttot 

(xlO"^) 


Ttot 

(xlO"^) 


Ttot 


n 

p° 


775.49 
±0.34 


149.4 
±1.0 


1.3 


100%) 


4.7 


4.6 


6.0 X 10"^ 




782.65 
±0.12 


8.49 
±0.08 


23.2 


(89%) 


7.2 


9.0 


8.9% 




1019.455 


4.26 


46.2 




29.7 


28.6 


12.7 X 10"^ 




±0.020 


±0.04 




(49%) 









* These branching ratios are at the pole mass. 



HELIOS/3 and NA60 at CERN SPS overcame these dif- 
ficulties and successfully measured the dimuon spectra 
all the way down to the pair-mass threshold of 2m^.^ 



Poisson-distributed, and that the like sign pairs are un- 
corr elated, the combinatorial background B can be ac- 
counted for by 



2. e^e -pair detection 

The spectrometer used in the discovery of Jj^jcc) 
in the p ± Be ^ e+e~X reaction (jAubert et al\ , Il974l ) 
contains the essence of e+e~ measurement, such as (i) 
excellent electron identification (hadron rejection), (ii) 
good momentum (pair-mass) resolution, and (iii) impor- 
tance of rejecting e+e~ pairs from photon conversion and 
Dalitz decays. 

In the e+e~ spectra, severe background sources are 
photons from meson decays such as tt^ 27 converting 
in the target (7 e+e~) and in detectors, and the Dalitz 
decays such as tt^ 7e+e~, 77 7e+e~, rj' ^e^e~ 
and uo 7r^e+e~. 

Although pairs from these sources have small opening 
angles and low masses, the limited track reconstruction 
efficiency and acceptance lead to a combinatorial back- 
ground for events in which two or more of these low-mass 
pairs are only partially reconstructed. This is the central 
problem of any low-mass e+e~-pair experiment. 



3. Combinatorial background 

Uncorrelated sources can produce, in addition to 
unlike-sign pairs, like-sign (/+/+ and /~/~) pairs. 

Most experiments make use of this fact in evaluating the 
combinatorial background. 

A typical method of subtracting the c ombinatorial 
background is as follows (jToial . l20od l2007h : Under the 
assumption that electron and positron multiplicities are 



Since muons must be energetic enough to penetrate the hadron 
absorber, the rapidity of the reconstr ucted pairs a r e high , e.g., 
3.3 < y < 4.3 in the case of NA60 (iDamianovid . 120071) . For 
e+e~, pa irs can be measured in t he mid-rapidity region of 2.1 < 
J] < 2.65 (lAgakichiev et aLLl2005h . 



B = 2VA/'++A/'— , (V.80) 

where A/'++ and N are the number of measured 
and pairs, respectively. The number of signal pairs 
S is then obtained as = A^^ — 5, where A^+~ is the 
number of measured unlike-sign pairs. 

This would work if the detector acceptance is the 
same for like and unlike sign pairs, and if a sufficient 
number of like-sign pairs are collected. This is in gen- 
eral not the case. A mixed event technique is then 
used to compute the combinatorial background. In this 
method, unlike-sign tracks from different events (with 
similar event topology) are paired. Since the same track 
can be used many times, paired up with tracks from dif- 
ferent events, the background spectra can be generated 
with high statistics. The accuracy of the technique can 
be tested by comparing the shape of the measured like- 
sign pair spectrum with that of the mixed combinatorial 
background. The generated background event distribu- 
tion can be normalized to the number of expected unlike- 
sign pairs from Eq. (|V.8Qp . 



C. High Energy Heavy Ion Reactions 

Although the main subject of the present review is the 
behavior of mesons produced in nuclei with elementary 
reactions, we nevertheless touch upon the low- mass (^ 
IGeV/c^) dileptons observed in heavy- ion collisions. 

An enhanced yield of dilepton pairs in the low-mass 
region is ubiquitous from 1 AGeV (Bevalac/SIS), through 
SPS energies (40 - 200 AGeV), up to the RHIC energy of 
^^s^ = 200 GeV. Here, the "enhancement" is defined 
as the excess of the observed yield over the sum of 
the "hadronic cocktail" as discussed above. 
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FIG. 36 Top view of the DLS (d ilepton spectrometer) at Be- 
valac (|Yegneswaran et oil . llQQOl ). 
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FIG. 37 The dilepton spectrum for Ca+ Ca at 1.0 AGeV mea - 
sured by the DLS collaboration (circles) ([Porter et all Il997l ) , 
compared with th e "hadronic cocktail' ' assum ing the "free" p 
spectral function ([Bratkovskava et a/.UlQQSh . 



1. Bevalac/SIS energies (1 - 2AGeV) 
DLS 

The first anomalous dilepton excess was reported by 
the DLS (dilepton spectrometer) experiment at Bevalac 
(lAGeY). Using a two-arm spectrometer as shown in 
Fig. [361 with a pair-mass resolution of AM/M ~ 10%, 
they succeeded for the first time to measure the dielec - 
tron spectra in heavy-ion collisions ([Porter et al 1 I1997D . 
When the measured spectra were compar ed with trans- 
port theory calculations ( Bratkovskava and Gassing , 



2008 : Gassing and Bratkovskavai 1 19991: iCozma et al 



20061 : lErnst et a/.LIl998l : IShekhter et a/.Ll2003h . an excess 
of about a factor 6 — 7 was found in the mass range of 
0.15 < M^^^_ < 0.4 GeV (Fig. [371). Including the p- 
meson modifications in the medium (mass dropping) did 
not eliminate the discrepancy (the observed yield was 
still higher by about a factor of 3 over the HSD curve). 
This has become known as the "DLS puzzle" . 




FIG. 38 A side view of HADES. The RICH detector, consist- 
ing of a gaseous radiator, a carbon fiber mirror and a tilted 
photon detector, is used for electron identification. Two sets 
of multiwire drift chambers (MDCs) are placed in front and 
behind the magnetic field to measure particle momenta. A 
time of fiight wall (TOF/TOFINO) accompanied by a Pre- 
shower detector at forward angles is used for an additional 
electron identification and trigger purposes. For a reaction 
time measurement, a start detector is located near the target 
([Salabura et a/.Ll2004l ). 
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FIG. 39 Direct comparison of the dilepton pair mass distribu- 
tions measured in C + C at 1 AGeV by HADE S (within the 
DLS acceptance) and at 1.04 AGeV by DLS ([Porter et all 
ll997l V Statistical and systematic errors are shown. Overall 
normalization errors (not show n) are 20% for the HADES and 
30% for the DLS data points ([Agakishiev et al 1 120081 ). 



HADES 

This DLS puzzle has been recently revisited by the 
HADES (high acceptance dielectron spectrometer) col- 
laboration at the heavy ion synchrotron SIS at GSI 
Darmstadt. HADES uses modern technologies such as a 
ring-imaging Cerenkov detector (RICH) to achieve good 
particle identification as well as a high mass resolution of 
AM/M 2.7% (Fig. [38]). They recently studied lAGeV 
C + C collisions in a low-resolution mode {AM/M = 8% 
at 0.8 GeV/c^ to emulate that of DLS) and projected the 
measured spectra into the DLS acceptance. 
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FIG. 40 Helios-3 setup (|Angelis et a/.L 120001 1. 
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The resulting spectrum, shown in Fig.[39l is consistent FIG. 41 M^^ distribution normalized to the charged-particle 



with that measured by DLS. For masses of 0.15GeV/c < 
M ^ _ < 0.50GeV/c^ it exceeds expectations based on 
the known production and decay rates of hadrons (most 
important being the rj meson) by a factor of about 7, 
thereby reconfirming the DLS data. 

However, recent HADES measurements of p -\- p 
e^e~X and p + n e^e~X seem to show that 
the C + C spectrum at lAGeV agree well with the 
\{pp + n p) spectrum at 1 .25 G eV, when scaled by the 
TT^ yield ([Galatvuk et al 1 12QQ9D . This may be indicat- 
ing that the DLS/HADES effect does not have nuclear 
(in-medium) origin, but that the NN bremsstrahlung 
cross s ections are in fact larger than hitherto as- 
sume d ( Kaptari and Kampferl . l2006l : IShvam and Mosel . 

With the enhanced bremsstrahlune; cross 
sections implemented in the HSD transport code, 
iBratkovskava and Cassin"3 (|2008l ) have recently shown 
that the calculated spectra agree well with DLS and 
HADES data. 



multiplicity obtained by the HELIOS/3 collaboration for S- 
W (fill ed squares) and p-Y^ (open squares) collisions at 200 
AGeV (jAngelis et all l2000l ) . 




FIG. 42 Schematic view of the CERES spectrometer with a 
radial drift time projection chamber (TPC). The latter was 
added in NA45/2 for the Pb+Au runs (,Marfn. . .2004 ) . 



2. SPS energies (40 - 200AGeV) 

At the CERN Super Proton Synchrotron (SPS), low- 
mass dilepton spectra were studied in the e+e~ mode 
by the CERES collaboration, and in the mode by 

the HELIOS/3 collaboration and the NA60 collabora- 
tion. These experiments all reported a low-mass dilepton 
enhancement. 



HELIOS/3 (m^m") 

The HELIOS/3 experiment used a dimuon spectrom- 
eter shown in Fig. [40] to measure distributions 
in proton on tungsten and sulphur on tungsten at 200 
AGeV . The spectrometer consisted of a hadron absorber 
(placed 25 cm downstream from the target), six interac- 
tion lengths of AI2O3 and 100 cm of Fe, followed by a 
magnetic spectrometer and muon hodoscopes. 

They found, by comparing the measured p-W and S- 
W dimuon distributions, each normalized to the charged- 
particle multiplicity (Fig. [41]), an excess in S-W interac- 
tions relative to minimum-bias p-W interactions. The ob- 



served excess is continuous over the explored mass range 
and has no apparent resonant structure. In the low mass 
(< 0.7 GeV) region the dimuon yield increases by 76±4% 
of the corresponding p-W dimuon spectrum (in the higher 
mass region, the excess was higher). 



CERES (NA45) (e+e") 

CERES is an innovative 'hadron-blind' axial- 
symmetric detector (Fig. [42]) dedicated to the measure- 
ment of electron pairs in the low-mass range (up to 
-1.5 GeV/c^). At the heart of CERES are the two 
coaxial ring-imaging Cerenkov detectors (shown in the 
left half of Fig. [42]) having a high Cerenkov threshold of 
7thr — 32, placed in a superconducting solenoid. With 
this setup, e+e~ pairs in a window of A77 = 0.53 around 
mid-rapidity were selectively detected and reconstructed. 

The dielectron pairs per charged particles, shown in 
Fig. [43] show that while the p-Be and p-Au data are 
reproduced within errors by Dalitz and direct decays 
of neutral mesons as known from p-p collisions, dielec- 
trons from S-Au collisions reveal a substantial enhance- 
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FIG. 43 CERES inclusive e+e" mass spec tra of 450 GeV p- 
Be, p - Au, a nd 200 AGeV S-Au collisions (|Agakichiev et all 
Il998l . Il995l ). Plotted is the number of electron pairs per 
charged particle, both in the acceptance and per event. Con- 
tributions from various hadron decays as expected from p-p 
collisions are shown together with their sum (thick line), and 
the s ystematic error on the latter is indicated by the shaded 
area ([Adamova et a/.Ll2008l V 
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FIG. 44 Layout of the NA 60 detectors in the vertex region 
(|Usai. G. and othersl . [20051 ) . 



m ent in the mass r egion 0.2 — 1.5GeV/c^ of a factor 
5 (jAgakichiev et a/.LIl995h . 

This observation generated lots of excitement in the 
community. It has been attributed to the pion anni- 



hilation in the fireball, tt+tt 



with 



strong in-medium mo dification of the interm ediate p, 
such as mass dro pping ("Brown an d Rhd l2002l ) or broad- 
ening ffiapp and Wambach. , ■200Q) . 



NA60 

The NA60 experiment added a telescope of radiation- 
tolerant silicon pixel detectors in between the target and 
the hadron absorber of the NA50 dimuon spectrometer 
(see Fig. H^ . This enabled the collaboration to match 
muon tracks before and after the hadron absorber, both 
in angular and momentum space, thereby improving the 
dimuon mass resolution in the re gion of light vector 
mesons fron i ^ 80 to - 20MeV/c2 ()Arnaldi et all l2QQ6l : 
lDamianoviAl2QQ7h . 

With this setup, NA60 succeeded to completely resolve 
uj and (p peaks in the In-In collisions at 158 AGeV^ for the 
first time in nuclear collisions. This is shown in Fig. 25] 
(the mixed-event technique was employed here to sub- 
tract the combinatorial background). 

By adjusting the cross section ratios r]/uj^ p/uo and 
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FIG. 45 Isolation of an excess above the hadron decay cock- 
tail (see text). Total data (open circles ), individual cock- 
tail sources (solid ), difference data (thick tria ngles ), sum of 
cocktail sources and diff erence data (dashed ) (|Arnaldi et all 
I2OO6I : lDamianovi3 . l2007l ) . 




FIG. 46 Comparison of the excess mass spectrum for 
the semi-central bin to model predictions, made for In-In 
at dNch = dri\ri^o = 140. Cocktail p (thin solid red 
line) , unmodifie d p (dashed-do t ted r e d line), in-mediun i 
broadening p (|Chanfray et all Il996l : iRapp et~oIl Il997l : 
iRapp and Wambachl . I2OO0I) (thick s olid blue li ne), in-medium 
moving p related to (|Brown and R ho. 1993, 12002I : [l i et all 
I1995I ) (dashed green line). The errors are purely statistical. 
The systematic error s of the continuum a re about 25%. From 
lArnaldi et all (|2006l ): iDamianovid (|2008l ). 



as well as the level of D meson pair decays, the periph- 
eral data could be fitted by the expected electromagnetic 
decays of the neutral mesons, i.e., the 2-body decays of 
the 7^, p, uj and (j) resonances and the Dalitz decays of the 
77, T]^ and LO. In the more central cases, a fit procedure is 
ruled out due to the existence of a strong excess with a 
priori unknown characteristics. 

The excess was therefore isolated by subtracting the 
measured decay cocktail, without the p, from the data, 
as shown in Fig. [45l The resultant distribution shows 
some non-trivial centrality dependence, but is largely 
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FIG. 47 (left) CERES unlike-sign pair yield (histogram) and 
combinatorial background (dashed curve). (right) Invari- 
ant e^e~ mass s pectrum compared to the expectation from 
hadronic decays (|Adamova et all l2008l ) . 



consistent with a dominant contribution from p 
fi^ li~ annihilation. Fig. [46] shows a distribution ob- 
tained for semi-centr al collisions, compar e d wit h in- 
medium broad ening (iRapp and WambachI , |2QQQ| ) and 
mass-dropping ferown and Rhol. l2QQ2h scenarios. The 
observed distribution {p spectral function) exhibits con- 
siderable broadening, but essentially no shift of the p- 
peak position. 



CERES (NA45/2) (e+e") 

In preparation for the lead beam acceleration in the 
SPS, CERES upgraded the detector by adding a cylin- 
drical time projection chamber (TPC) with a radial elec- 
tric field (right half of Fig. H2|) . Among other things, 
this improved the mass resoluti on Am/m in the region 
of the p /cj from 9% to about 6% (jAgakichiev et a/.L[2QQ5l : 
iMarinl . [200 



I2QQ4 ). 



The dielectron distribution obtained in Pb-Au colli- 
sions at 158 AGeV before combinatorial subtraction, to- 
gether with the normalized mixed-event background is 
shown in the left panel of Fig. SJl The background- 
subtracted distribution is compared with the hadronic 
cocktail in the right panel of Fig. [371 Here again, an 
enhancement over the cocktail is observed in the mass 
range 0.2 < m ^ _ < l.lGeV/c^, the enhancement fac- 
tor being 2.45 ±0.21 (stat) ±0.35 (syst) ± 0.45 (decays), 
where the last error is from the systematic uncertainty 
in the cocktail calculation. 

Fig. [IHl shows the dielectron yield after the hadronic- 
cocktail subtraction, compared with the mass dropping 
(left) and width broadening (right) assumptions. Al- 
though the error bars are larger than those in the NA60 
spectra, the authors concluded that a substantial in- 
medium broadening of the p is favored over a density- 
dependent shift of the p pole mass. 



FIG. 48 CERES e^e~ pair yield after subtraction of the 
hadronic cocktail. In addition to the statistical error bars, 
systematic errors of the data (horizontal ticks) and the sys- 
tematic uncertainty of the subtracted cocktail (shaded boxes) 
are indicated. The broad e ning scenario (long-d a shed line: 
(Ivan Hees and Rapd. l2006l : iRapp and Wa mbachl . l2000l )) is 
compared to a calculation assum ing a density dependent drop- 
ping p mass (dotted line in (a): ([Brown and Rhd. 119911 . [19961 . 
[2003)) and to a broadenin g scenario excluding b aryon effects 
(dotted line in (b)). From lAdamova et al\ (|2008l ). 




FIG 49 The PHENIX experiment at RHIC (jAdcox et all 
I2OO3I ). 



3. RHIC 



200 GeV) 



At RHIC, the PHENIX experiment (Fig. [49]) has been 
designed to measure dielectrons over a wide mass range. 
Electrons and positrons are reconstructed in the two cen- 
tral arm spectrometers using drift chambers, located out- 
side an axial magnetic field. They are identified by hits 
in the ring imaging Cerenkov detector (RICH) and by 
matching the momentum with the energy measured in 
an electromagnetic calorimeter. 

Fig. [50] shows a dielectron distribution observed by 
PHENIX in Au+Au minimum-bias collisions at ^/s^ = 
200GeV, after combinatorial background subtraction us- 
ing the mixed event shape norm alized to the like-sign pair 
yields ([Afanasiev et all [2007h . The dielectron yield in 
the minimum bias collisions, in the mass range between 
150 and 750 MeV/c^, is enhanced over the cocktail by 
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FIG. 50 Invariant e^e~ pair yield observed by PHENIX in 
Au+Au minimum-bias and in p-\-p collisions at ^s^^ = 200 
GeV. The curve is the hadronic cocktail for p -\- p. The 
p + p spectrum and the cocktail curve were normalized 
to the Au+Au yield in the ^ + < 100 MeV/c^ region 
(jAdare et all l2008l : lAfanasiev etall l2007l : iToial . 120081 ). 



a factor of 3.4 ± 0.2(stat.) ± 1.3(syst.) ± 0.7(model). A 
clear increase with centrality is also observed. No de- 
tailed analysis of the excess is available yet. 



4. High energy heavy ion summary 

An enhanced yield of dilepton pairs over the hadronic 
sources in the low-mass region has been observed, regard- 
less of the bombarding energy. 

At low beam energy of ~ lAGeV, the long-standing 
DLS puzzle (excess) has been confirmed by the recent 
HADES experiment. However, recent indications of C+C 
dielectron distribution agreeing with the {pp-\-pn)/2 dis- 
tribution, if confirmed, may rule out the possibility of 
in-medium modification effects at this energy. 

The two SPS experiments, CERES measuring dilep- 
tons and NA60 measuring dimuons, both established that 
there is a dilepton enhancement in the low-mass region. 
The excess here is consistent with tt+tt" p 
with the p significantly broadened in the nuclear medium, 
while the data do not call for the simple p mass change. 

At RHIC, the PHENIX experiment showed that there 
is a dilepton enhancement in the low-mass region, the 
magnitude of which increases faster with the centrality of 
the collisions than the number of participating nucleons, 
but the statistical errors are still fairly large in order to 
draw firm conclusions based on the data. 




FIG. 51 A pla n view of the TAGX spectrometer 
(|Huber ali l2003l ). 



University) used a tagged photon beam in the energy 
range of 600-1120 MeV to study the 7A tt+tt-X re- 
action on A =^H, ^He and ^^C targets. This was a pio- 
neering experiment which attempted to study in-medium 
modifications of p^ with elementary reactions. 



However, the claim of 
-160 ± 35MeV/c2 in ^He 



finding 



a 



mas s shift of 
19981 ) was met 



D. p, uj and 
reactions 



mesons produced in nuclei with elementary 



m "JrLe (|Lolos et al 
with skepticism due to the inevitable pion rescattering 
effect even for light targets (the emitted pions being in 
the resonance region), the small target volume, and the 
much-larger-than-expected shift. 

They later applied a helicity analysis t o extract in- 
medi um invariant mass distributions ( Huber et all 
l20Q3h, conipared the spectra with various simulation s 
(jPost et all I2OO4I : iRaPD et all Il997l : ISaito et a/.LIl997af ). 
and obtained a smaller but still sizable mass shift of 
—65 ~ — 75MeV/c^ in the photon energy bin of = 
800 - 900 MeV, and -45MeV/c2 for = 960 - 
1120MeV. The ^^C distributions on the other hand, were 
found to be consistent with quasi-free p^ production. 

Why TAGX observed such a large effect only in ^He is 
not yet understood, but in view of the fact that the 
e~^e~X data of CLAS-g7 (section rV.D.3p do not show any 
sign of p^ mass shift, this is most likely unrelated to the 
p^ in- medium modification.^ 



2. E325 experiment at KEK 

The experiment E325 at the KEK 12 GeV Proton Syn- 
chrotron was the first to measure dileptons in search for 
the modication of the vector meson mass in a nucleus in 
elementary reactions. They measured the invariant mass 
spectra of e+e~ pairs produced in 12 GeV proton-induced 
nuclear reactions. The setup is a two-arm spectrometer 
(Fig. [52]), and was designed to measure the decays of 



1. TAGX at INS Electron Synchrotron 

The TAGX experiment (Fig. [5lD at the 1.3-GeV INS 
Electron Synchrotron (Institute for Nuclear Study, Tokyo 



^ The incident photon energy range of TAGX was 0.6 — 1.12 GeV 
while it was 0.61-3.82 GeV in the CLAS-g7 experiment. There- 
fore, a more direct comparison would be to use low incident en- 
ergy events of the CLAS-g7 data sample. 
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FIG. 52 Schematic view of the experimental setup of the E325 
spectrometer: (a) the top view and (b) the side view. The 
side view shows the cross section along; the cent er of the kaon 
arm (|Ozawa et ali l2QQll : ISekimoto etal\ , 120041 ) . 




2 3 

rapidity 



FIG. 53 Kinematical distributions of e^e pairs in the cj) mass 
region (0.95 < 1.05GeV/c^ detected in the E325 

spectrometer (points with error bars) , together with the simu- 
lation result using the JAM nuclear cascade code (histogram, 
(|Nara et a/.L 119991 )). (left) p-f distri bution, (right) R apidity 
y vs pair transverse momentum pr (|Muto et all l2007l ) . 



TABLE II Expected in-nucleus decay fractions of vector 
mesons in the E325 kinematics, assuming that the meson 
decay widths are unmodified in nuclei, obtained by 
a Monte Ca r lo-typ e model calculation 
iNaruki et all . l2006f ). 



obtained by using 
(iMuto et ali l2007l : 





c 


Cu 


p 


46% 


61% 


ix) 


5% 


9% 







6%' 



for slow (j) mesons with < 1.25. 



the vector mesons, ^ e+e , p/cj ^ e+e as well as 
^ ^ K+K-. 

For electron identification, two stages of electron- 
identification counters were used. The first was the front 
gas-Cerenkov counters (FGC). The second stage con- 
sisted of the rear gas-Cerenkov counters (RGC), the rear 
lead-glass electromagnetic (EM) calorimeters (RLG), the 
forward lead-glass EM calorimeters (FLG), and the side 
lead-glass EM calorimeters (SLG). The overall electron 
efficiency was 78% with a pion rejection power of 3 x 10~^ 
(jSekimoto et a/.l . l2004D . 

The mass resolution was estimated to be 8.0 MeV/c^ 



and 



e+e decays. 



and 10.7 MeV/c^ for co 
respectively. 

The kinematical region covered was 0.5 < ^ < 2 and 
1 < /37 < 3 for e+e~ pairs (Fig. [53]), where the decay 
probability inside the target nucleus was expected to be 
enhanced. Assuming that the meson decay widths are 
unmodified in nuclei, the coverage would correspond to 
the in-nucleus decay fractions shown in Table [Til 

The E325 invariant mass spectra for C and Cu tar- 
gets are shown in Fig. [SH The data were taken with the 
"unlike-sign-double-arm" trigger condition, i.e., either a 
positron in the left arm and an electron in the right arm 
(LR event) or vice versa (RL event), to suppress the back- 
ground from Dalitz decays and conversions, thereby pre- 
cluding the possibility of normalizing the background to 
the like-sign pair distribution. 



The combinatorial background shape was obtained by 
the event-mixing method, and its normalization was ob- 
tained by fitting the data together with contributions 
from CO e+e~, p e+e~, (j) e+e~, r] e+e~7 
and CO e^e~7r^. The relativistic Breit-Wigner distri- 
bution was used for the resonance shapes, and kinemati- 
cal distributions of mesons were obtained by the nuclear 
cascade code JAM ( Nara et all Il999l ). which is in good 
agreement with the experimental data (see, e.g.. Fig. [53] 
(a)). 



E325 results on the p/co mesons 

The striking features of the p/co region of the E325 
spectra are as follows, i) A significant excess can be 
seen on the low-mass side of the co peak, which could 
not be fitted with the cocktail. Therefore, in the fit 
shown in Fig. [54] the range 0.6 < m ^ _ < 0.76GeV/c^ 
was excluded from the fit. ii) The p/co ratio, which is 
known to be close t o unity in pp collisions at this energy 
(|Blobel all 11974 ). is here p/co < 0.15 and < 0.31 for 
C and Cu targets, respectively (95% C.L.). 

The disappearance of p and the appearance of the ex- 
cess may be due to the in-medium dropping of the p 
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FIG. 54 Invariant mass spectra of e^e~ for the (a) C and (b) 
Cu targets. The sohd hnes are the best-fit results, which is 
the sum of the known hadronic decays, oo e^e~ (dashed 
line), (j) e^e~ (thick dash-dotted line), 77 e^e~7 (dash- 
dotted line), and to e^e~7T^ (dotted line) together with the 
combinator ial background (lon g- dashed line), p e^e~ is 
not visible (iNaruki a/.Ll2006l ). 
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FIG. 55 Invariant mass spectra of 

background and the shapes of 77 e^e~7 and 00 ^ e' e tt 
were subtracted. The result of the model calculation consid- 
ering the in-medium modification for the (a) C and (b) Cu 
targets. The solid lines show the best-fit results. In (a) and 
(b), the shapes of a; ^ e^e~ (dotted line) and p e^e~ 
(dash-dotted line) were modified according to the model us- 
i ng the formula my ip) Imyj^) = 1 — k{p/po) with k — 0.092 
(jNaruki et a/.l . l2006l ) . 



mass^, which would take away the strength from the nor- 
mal p and put them in the excess region. This assumption 



^ Alternatively, this may be due to the over-subtraction of the 
background component, as pointed out by the J-Lab CLAS g7 
collaboration, whose results are in conflict with those of KEK 
E325. See Section EES 



was tested by fitting the background-subtracted spec- 
tra using a Monte Carlo-type model including the mass- 
dropping model, 



my(p)/my(0) = 1 - k{p/po). 



(V.81) 



The vector mesons were generated on the surface of 
an incident hemisphere of the target (supported by 
the A^^^ dependence of the uj production cross section 
( Tabaru et al propagated through the nucleus 

which was modeled by a Woods-Saxon density distribu- 
tion. The parameter k was common for co and p as well 
as for C and Cu targets. The p/uo ratio was also allowed 
to vary. 

The best fit results are, k = 0.092 ± 0.002 and a p/o; 
ratio of 0.7 ± 0.1 and 0.9 ± 0.2 respectively, for C and 
Cu targets. The best-fit curves are superimposed on the 
background-subtracted spectra shown in Fig. [55l They 
also examined whether or not the p — uo interference can 
account for the observed shoulder, but found that the in- 
terference cannot explain the data even though the p/uo 
ratio and the mixing angle were scanned over a wide 
range. 

So, for the p/uo region, E325 concluded that both p and 
uo masses are shifted by 9% at the normal nuclear density. 
Fits with density-proportional width broadening did not 
fit th e data; the fit resul ts favored the zero-broadening 
case ( Naruki et al ], l2006[ ). This is in conflict with the 
J-Lab CLAS g7 result discussed in lV.D.3[ 



E325 results on the meson 

Due to the long lifetime of the (j) meson, in-medium 
modification effects, if any, are expected only for the slow- 
moving mesons which have a chance to decay inside the 
target nucleus. E325 therefore divided the data in three 
parts based on the jS^ values of the observed e+e~ pairs, 
/37 < L25, L25 < p-f < L75 and L75 < (see Fig. [53]). 

The /37-selected spectra are shown in Fig. [56l together 
with fit results. The (j) was assumed to have in- vacuum 
mass and width, convoluted over the detector response in 
the simulation according to the JAM-generated kinemat- 
ical distributions of the (f) meson in each [5^ region. A 
quadratic background was added to the simulated peak, 
and the background parameters and the (j) abundance 
were obtained from the fit. 

The fits are satisfactory, except for the (5^ < 
1.25 region of Cu data, in which a large excess of 
A/gxcess/ (A^excess + A^</)) = 22% was fouud. If this excess is 
to be ascribed to the in-medium (f) modification, not only 
the mass but also the width need to be varied, since the 
JAM-based simulation indicates only 6% of the (j) meson 
produced in copper nuclei with < 1.25 would decay 
in the target nucleus if broadening is not introduced (see 
Tab. ini). 

It was thus attempted to fit the data by introducing 
both the density-linear mass shift 

m^{p)/m^{0) = 1 - ki{p/po) 
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FIG. 56 Obtained e^e~ distributions with the fit results. The 
target and [3^ region are shown in each panel. The points 
with error bars represent the data. The solid lines repre- 
sent the fit results with an expected cj) e^e~ shape and a 
quadrat ic background. Th e dashed lines represent the back- 
ground (|Muto et a/.L 120071 ). 
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FIG. 57 Confidence ellipsoids for the modification parameters 
ki and k^^^ in cases (i) in (a) and (ii) in (b). The values 
of Ax^s in both panels are the differences from the Xmin(— 
316.4) at the best-fit point in case (i) which is shown by the 
cross in the panel (a). The best-fit point in case (ii) is shown 
by the closed circle in the panel (b), and also in (a) since 
the ordinate s are common to both cases in parameter space 
(|Muto a/.L 12007. ). 

and the density-linear width broadenings 

Tf\p)/rf\0) = l + krip/po), and 
ri%p)/Ti%o) = 1 + krip/po), 

where r^°* is the total width and is the e+e~ partial- 
decay width. 

Figure [53 shows the fit result for the two cases ex- 
amined, (a) = /c2°* (i.e., the branching ratio F^^/F^^* 
remains unchanged in the medium), and (b) /cf^ = (i.e., 
ptot increases but F^^ does not increase in the medium). 

The fit favors the former case. The obtained values 
are ki = 0.034^0.007 and /cf * = fcf = 2.6tli, indicating 
the in- medium 0- meson mass shift of 3.4% and width 




FIG. 58 (top) Side view of the CLAS detector in Hall B of 
J-Lab, with the photon-tagging s ystem, (bottom) CLAS cut- 
away view (Mecking et aa 120031 ). 

mcrease of a factor of 3.6 {Tf^ ISMeV/c^) at normal 
nuclear density. 



3. J-Lab EOl-112 (g7) experiment 

J-Lab EOl-112, better known as the CLAS experiment 
g7, was conducted in in Hall-B of Jefferson Laboratory 
(Fig. [58]). An electron beam accelerated by the Contin- 
uous Electron Beam Accelerator Facility (CEBAF) was 
used to produce a tagged photon beam having an energy 
range of 0.61 - 3.82GeV. 

The CEBAF Large Acceptance Spectrometer (CLAS) 
is a nearly 47r-detector based on a six-coil superconduct- 
ing toroidal magnet, and was designed to track charged 
particles with momenta greater than 200 MeV/c over the 
polar angle range from 8° to 142°, while covering up 
to 80% of the azimuth. The CLAS detector is divided 
into six identical spectrometers (sectors), each made of 
three regions of drift chambers (DC), time-of- flight scin- 
tillators, Cerenkov counters (CC) and electromagnetic 
calorimeters (EC) (Fig. [58] bottom). The target materi- 
als were liquid deuterium (LD2), carbon, titanium, iron. 
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TABLE III Mass and width of the p meson obtained by the 
CLAS-g7 collaboration from the simultaneous fits to the mass 
spectra for each target and the ratio to ^H. 



300 



200 



Target 




100 




773.0 ±3.2 185.2 ±8.6 




C 


726.5 ±3.7 176.4 ±9.5 





Fe, Ti 


779.0 ±5.7 217.7 ± 14.5 





and lead (simultaneously in the beam). To reduce the 
low-energy e~ and e+ background from pair production 
in the targets, a "mini-torus" magnet was situated just 
beyond the target region and inside the DC. 

The e+e~ event selection and the rejection of the very 
large tt+tt" background were done through cuts on the 
EC and the CC. The pion rejection factor was 5.4 x 10~^ 
per track, or 2.8 x 10~^ for the pair. The pair-mass reso- 
lution was 10 MeV/c^ for the cj) peak. The pair momenta 
of ~ 0.8 ~ 1.8Ge V/c were acce pted, similar to the KEK 
E325 acceptance (jPialalil . [2008h . 

In reconstructing the e+e~ pairs, the two leptons were 
required to be detected in different sectors of the CLAS 
detector. This requirement removed the large back- 
ground due to pair-production, Bethe-Heitler processes, 
and TT^ and r] Dalitz decays that have a small opening 
angle. 

The combinatorial background^ was approximated 
by an event-mixing technique, and was normalized to 
the number of expected opposite-charge pairs, calcu- 
lated from the number of observed like-sign pairs using 
Eq. (|V.80p . The spectra shown in Fig. [59] (left) are the re- 
constructed e+e~ distributions, co mpared with the nor - 
malized combinatorial background ( Wood et alV 120081 ). 
The uncertainty of the normalization was estimated at 
±7%. 

Monte Carlo calculations using a code based on a 
semi-classical BUU transport model were used to fit 
the background-subtracted spectra (Fig. [59l right). In 
the model, the particles produced as a result of the 
7A/' reaction in the target nucleus were propagated 
through the nucleu s allowing for final-state interactions 
( Effenberger et all , Il999l ). The acceptance-corrected 
BUU mass shapes for the p, uo and (j) mesons were scaled 
separately to match the experimental mass spectra. A 
substantial contribution from the p meson was found 
(dot-dashed curves in Fig. [59l right) unlike in the KEK 
E325 analysis. 

As the probabilities of the uo and (j) mesons decaying 
inside the nucleus are low, the simulated uj and (j) mass 
shapes were subtracted from the data, to obtain the p 
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FIG. 59 (left) Normalized combinatorial background for 
individual targets compared with the data. (right) Re- 
sult of the fits to the e^e~ invariant mass spectrum ob- 
tained fo^_the_^H (top) , C (middle), and Fe-Ti (bottom) 
data (|Wood et a/.L l2008l ). The curves on the right panels 
are Monte-Carlo calculations by t he BUU model fo r vari- 
ous vector meson decay channels (|Effenberger et al ], I1999I : 
■Effenberger and Mosel . 2QQQ1 



mass spectra. The results are shown in Fig. [6OI The 
curves therein are Breit-Wigner//i^ fits (/i being the in- 
variant mass). Here, the /i^ factor comes from the mass 
dependence of the r^^^_ (/i)/rtot(M) ratio; r^^^_ (/i) ex 
l/Zi"^ and rtot(M) pi. 

The mass and width of the p meson in various targets 
were obtained by performing a simultaneous fit to the 
mass spectra and the ratio of each spectrum to the 
data (so as to impose more constraints on the fits), and 
the fit results are shown in Table Hill These are consistent 
with collisional broadening without mass modification. 

The mass shift coefficient k as defined in Eq. ()V.8ip 
was obtained by analyzing the ratio of the Fe-Ti to the 
distributions to be 0.02 ± 0.02, which corresponds to 
an upper limit oik = 0.053 with a 95% confidence level. 
These results are quite different from those obtained by 
the KEK E325 experiment. 



^ In the CLAS g7 experiment, the probability of an untagged pho- 
ton and a tagged photon being in the same radio-frequency tim- 
ing bunch was about 25%. This contributed to the combina- 
torics, in addition to the usual case of picking up a wrong lepton 
produced in the same event. 



4. CBELSA/TAPS experiment 

The CBELSA/TAPS collaboration at the electron 
stretcher accelerator (ELSA) in Bonn used the 
n^^X reaction to study the uo meson in-medium behav- 
ior using the Crystal Barrel (CB) and TAPS crystal spec- 
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FIG. 60 Individual Breit- Wigner //x^ fits to the p mass spectra 
(background and cj, (j) contributions subtracted) (|Wood et a/.L 

mm- 
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FIG. 61 Side view of the Crystal Barrel (CB) and TAPS 
detector combination. 



trometers shown in Fig. [611 Tagged photons in the energy 
range of 0.64 — 2.53 GeV were incident on targets (Nb and 
LH2) mounted in the center of the CB, a photon calorime- 
ter consisting of 1290 CsI(Tl) crystals with an angular 
coverage of 30° to 168° in the polar angle and a complete 
azimuthal angle coverage. Reaction products emitted in 
forward direction were detected in the TAPS detector, 
which consisted of 528 hexagonally shaped BaF2 detec- 
tors covering polar angles between 4° and 30° and the 
complete 27r azimuthal angle. The resulting geometrical 
solid angle coverage of the combined system was 99% of 
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FIG. 62 (left) The mass distribution obtained from a 
Monte Carlo simulation of the process 7 + Nb 7r°7 + X 
at = 1.2 GeV. The spectrum is decomposed into differ- 
ent contributions corresponding to the fraction of cj-mesons 
decaying outside the nucleus (a), the fraction of cj-mesons de- 
caying inside for which the 7r° does not rescatter (b), and the 
fraction of cj-mesons decaying inside the nucleus for which 7r° 
rescatters(c). In the simulation a drop of the uo mass by 16% 
at normal nuclear density was assumed, (right) The same 
as the left panel, with the additi onal condition of T^o > 150 
MeV (jMesschendorp et ali l200lh . 



47r. Charged particles were identified with a scintillating 
fiber detector placed inside the CB, and a pl astic scintil- 
lator mounted in front of each TAPS cryst al (|Aker et a/.L 
I1992I : IJanssen et all . I2OOOI : iNovotnvl . fT99lh . 

The uo 7r^7 decay mode has a large branching ra- 
tio of 8.9% and is a clean and exclusive mode to study 
the uo in- medium properties since the p ^ tt^j branch- 
ing ratio is only 6.0 x 10~^ (see Table |T)). Therefore, 
the study o f this mode is c omple mentary to the dilep- 
ton decays ( Sibirtsev et all I2OOOI ). A serious disadvan- 
tage are possible strong final-stage interactions of the 
TT^ meson within the nucleus. Monte Carlo simulations 



( Messchendorp et~al\ . l200l[ ) show that the rescatter ing 
effect is small in the mass range of interest, and can be 
further reduced by removing low-energy pions (T^^ < 150 
MeV), as depicted in Fig. [621 

The 7r^7 events are reconstructed from three photons, 
and the invariant mass spectra are shown in the left 
panel of Fig. [63l Here, in order to maximize the in- 
nucleus decay probability, slow-moving uo mesons with 
|pa;| < 0.5GeV/c were selected. The large continuum 
background is due to four-photon decays of tt^tt^ and tt^t] 
where one of the four photons is missed. A smooth poly- 
nomial background was assumed and was subtracted, and 
the resultant LH2 and Nb data are compared in the right 
panel. As shown, a shoulder on the low-mass side of the 
uo peak was found on the Nb target. This was taken as ev- 
idence for an uo in-medium mass reduction by 6OI35 MeV 
at an average nuclear density of O.Qpo^ or in terms of the 
mass shift coefficient k as defined in Eq. ()V.8ip . this gives 
k ^ 0.14. The width was found to be F = 55MeV/c^ 
dominated by the experimental resolution. 

The background- s ubtrac tion procedure was criticized 
by iKaskulov et all ( 20071 ), who pointed out if the 
same background shape is used both for LH2 and 
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FIG. 63 (a) Inclusive invariant mass spectra for momenta 
less than 500 MeV/c. Upper histogram: Nb data, lower his- 
togram: LH2 target reference measurement. The dashed lines 
indicate fits for the respective background, (b) 7r°7 invari- 
ant mass for the Nb data (solid histogram) and LH2 data 
(dashed histogram) after background subtraction. The error 
bars show statistical uncertainties only. The soli d curve repre- 
sents the simulated line shape for the LH2 target (|Trnka et all 
H). 
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FIG. 64 Mass number (A) dependence of the transparency 
ratio T. The dotted line corresponds to a a oc The 
solid and dashed cu rves show the theoretical calculation 
(|Cabrera et a/.L [20041 1 without (soUd curve) and with (dashed 
curve) Pauli-blocking correction for the meson scattering 
angle in the laboratory frame of 0° (;Ishikawa et a/.. , .2005 ) . 



Nb, the shoulder structure would disappear. The 
CBELS A/TAPS group pointed out that the experimen- 
tal data clearly show that the background distributions 
are different and hence it is not justified to assume the 
same background shape ( Metad . l2QQ8a| ). However, the 
fact that slightly different background assumptions lead 
to a complete different conclusion on the co mass shift is 
quite alarming. 

The CBELSA/TAPS group has therefore started to 
employ the event-mixing technique to generate the back- 
ground distribution, instead of using a polynomi al func- 
tion. Preliminary results were presented in (|Metagl , 
but these were later found to con tain some prob - 
lems, and are being further investigated (lMetad . [2QQ8bl ). 
Therefore, until the reanalysi s is finalized by th e group, 
the u mass shift reported in ( Trnka et all l2QQ5l ) cannot 
be regarded as a conclusive evidence for the in-medium 
CO modification. 



E. Vector-meson in-medium width from transparency-ratio 
measurements 

Instead of obtaining the in-medium meson width from 
fits to the observed invariant- mass peak, an alternative 
method of using the transparency ratio T, defined as 



for the loss of vector-meson fiux via inelastic processes 
in nuclei, and is related to the absorptive part of the 
meson-nucleus potential. 

This is conceptually a simple measurement, but ex- 
tracting the in-medium meson width from the A depen- 
dence of the ratio T requires comparison with theory cal- 
culations. 



The (j) attenuation 

The photo-produciton of (j) mesons from Li, C, Al and 
Cu targets was measured at = 1.5 — 2.4 GeV, using 
the laser-electron photon faci lity at SPring-8 (LEP S), in 
the jA K+K-X channel (jlshikawa et a/.Ll2QQ5[ ). The 
A dependence of the incoherent (j) photo production cross 
section was found to be cfa oc A°-72±o.07 ^ A~^-^^, 

as shown in Fig. [6lj). Using a Glauber- type model cal- 
culation, the in-medium ^-nucleon cross section was de- 
duced to be (j0Ar = 35^]^^ mb, which is much larger than 
the free-space value of a^^ff = 14Q/ib used as an input to 
the model calculation (* Amsler et all [2QQ8I ) . Theoretical 
calculations (.Cabrera et al.^ 200j) predicted much larger 
T values (solid and dashed curves in Fig. [64 

Using the classical low-density relation 



T = 



O'-fA^VX 



Act, 



'-fN^VX 



(V.82) 



was propos ed ([Hernandez and Oset', '19921; 

[Kaskulov etal] , [200/1 : Muhlic h and Mosel, 2006j, 
and ha s been used to extra ct (j) (jlshikawa et all [20051 ) 
and uj ( Kotulla et al I [2008[ ) in- medium widths. Here, 
a^A^vx is the inclusive nuclear vector-meson (V) 
photo-production cross section and ct^n^vx is the cross 
section on a free nucleon. The ratio T is a measure 



Fy = hpPca^ 



(V.83) 



this would correspond to a width of F^ :^ 80 MeV/c^ at 
p = po diiid p 0.7 (i.e., 1), where the KEK E325 

experiment reported a much smaller in-medium cj) width 
of 15 MeV/c^ (see section [EE!]) • 



This discrepancy may at least partly be due to the way the trans- 
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FIG. 65 Experimentally determined transparency ratio nor- 
malized to the carbon da ta in comparison with a theoretical 
Monte Carlo simu lation ([Kaskulov et all l2007l l (left) and a 
BUU calculation (iMuehlich et all 120041 ) (right) varying the 
width at 1.1 GeV/c momentum, respectively. The width is 
given i n the nuclear rest frame. Only statistical errors are 
shown (|Kotulla et all [2008|) • 



density p of the medium, and its experimental detection 
is one of the main challenges in modern hadron physics. 



Pion 

Theoretically, all hadrons receive various spectral 
changes due to their strong interaction with the medium. 
Among those, the in-medium modification of the pion de- 
cay constant /^(p) and /^(T) is theoretically well under 
contol at low T and p, and have close relation to the 
in- medium change of the chiral condensate pL42|IL43p . 
The predicted reduction at low density, (qq) ^ / {qq)^ — 
1 — (0.3 ~ 0.4)p/ po, has now been experimentally demon- 
strated by comparing the isovector pion-nucleon bi and 
pion- nucleus bi{p) scattering lengths derived from pio- 
nic hydrogen and deeply-bound pionic Sn atoms, respec- 
tively, as well as by analyzing the differential cross sec- 
tions for low-energy tt^ elastic scattering by several nu- 
clei. 



The LJ attenuation 



cr meson 



The CBELS/TAPS collaboration measured the A de- 
pendence of the CO photoproduction cross section on the 
nuclei C, Ca, Nb and Pb. The average momenta of the 
mesons was 1.1 GeV/c, so that almost all u mesons de- 
cay outside the nuclear target. Since the uj photoproduc- 
tion cross section on the neutron is not known, they took 
the transparency ratio normalized to the carbon data, as 
shown in Fig. [65l^^ 

The data were then compared with three different 
types of models, i) a Glauber model similar to the L EPS 
analysis, ii) a BUU analysis (|Muehlich et ah . 2004) and 
iii) a calculation by the Valencia group ([Kaskulov et all 
120071) . In all cases, the inelastic uj width was found to be 
130 — 150 MeV/c^ at p = po for an average uo momentum 
of 1.1 GeV/c, or in terms of uoN cross section, cf^^n — 70 
mb. 



VI. CONCLUDING REMARKS 



One of the interesting signals associated with the in- 
medium chiral restoration would be the spectral enhance- 
ment on the (J channel near the 27r threshold. Intriguing 
experimental results of "softening" of the (7r7r)/^j^o dis- 
tribution (i.e., shift of the peak position to lower masses) 
have been obtained. These agree fairly well with (i) in- 
medium modifications of the tttt interaction, as well as 
with (ii) rescattering of outgoing pions with the nucleons 

without in- medium t ttt interacti on. In fact, the results 

of two calculations, iRoca et a/] (|2002h and IBuss et al\ 
(|2006l ). predict very similar spectra, as shown in Fig. [35l 
However, as long as the rescattering scenario can repro- 
duce most of the observed "softening" trend, we cannot 
yet extract the predicted partial chiral restoration sig- 
nature from the tttt spectra. New high-statistics data 
taken on C, Ca, Pb with the Crystal Ball/TAPS detec- 
tor, which is beiM analyzed, may help shed some light 



on this problem (jMetag . 2008h ). 



The QCD vacuum shows the dynamical breaking of 
chiral symmetry. In the hot/dense QCD medium, the 
chiral order parameter such as {qq) (chiral condensate) is 
expected to change as a function of temperature T and 



parency ratio was normalized to the production cross section on 
the nucleon in Fig. I64|) . While the 77? cross section is sufficiently 
w ell known, the cross section is not. This was partly avoided 
in llshikawa et all (2005) by taking '^Li as a reference. 
For this reason, the ratio shown in Fig. [65] is named instead 
of T. This leads to different values of the transparency ratio since 
the transparency loss in C is normalized away. This normaliza- 
tion takes into account uj production processes on two nucleons 
which may be relevant in nuclei. 



Vector mesons 

Significant experimental work has been done to detect 
the possible in-medium "mass shift" of vector mesons, 
both using heavy-ion collisions and using elementary- 
particle beams. In general, the vector spectral function 
receives a shift of the peak, broadening, new structures, 
etc., due to the complex interaction of the vector current 
with the medium. Also, such a spectral shift may well de- 
pend on the spatial momentum of the current. Therefore, 
it would not be appropriate to oversimplify the problem 
to "mass shift vs. width broadening" . With this caution 
in mind, we list experimental results on the in-medium 
mass and width of the p, uo and (j) mesons produced with 
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FIG. 66 Comparison of the dielectron invariant- mass spec- 
trum (carbon target) of E325 (circles) and CLAS g7 (trian- 
gles). 



elementary reactions, measured in different experiments 
in Table IVl Tfie TAGX results are not included fiere for 
tfie reasons discussed in section IV.D.ll Tfie uo mass sfiift 
from CBELS A/TAPS is listed in tfie table, but it may 
change after the ongoing reanalysis, and hence we do not 
include this in the summary discussion. 

Upon examining this table, we realize that there are 
some inconsistencies, and we discuss the two most press- 
ing issues below. 

(i) E325 and g7 disagree on the p{u) mass shift: The 
E325 result is both p and uo masses get reduced at po 
by 9% (the mass shift parameter k = 0.092 ± 0.002), 
while the CLAS g7 placed a 95% confidence upper limit 
at /c = 0.053. The comparison of the background- 
subtracted dielectron distributions (carbon target) mea- 
sured by the two experiments (Fig. [66|) shows that the 
t wo spectra ar e very different 

IWood et all (|2008l ) pointed out that this difference 
must be due to the way the combinatorial background 
was subtracted in E325. In the E325 analysis, due 
to the lack of a sample of same-charged leptons by 
which to extract the normalization of the combinatorial 
background, the background contribution was fit along 
with the uj- and 0-meson shapes. Without an absolute 
determination of the combinatorial background, the 
p-meson signal was suppressed and included in the back- 
ground shape. Indeed, if the background normalization 
was free to vary in the CLAS-g7 fit, t he g7 spectra were 
found to be consistent with k ^ 0.16 (iDialalil . [20061 )^^. 



(ii) E325 and CBELS A/TAPS disagree on the co 
width: While an co width broadening was not observed 
by the E325 experiment, CBELSA/TAPS found an 



Note that the p/oo r atio in the pp c ollisions at the KEK energy 
is about unity ( Blobel et all 11974 ), while t hat in the 7j o colli- 
sions at the CLAS energy i s about 3 to 1 (iBarth et all l2003l : 
lErbe et all Il968l : IWu et all 120051 ). which may account for the 
bulk of the difference found in Fig. [66l 

In reality, there is no such freedom in the g7 background normal- 
ization. This was done just for the sake of g7-E325 comparison. 



unexpectedly large in-medium broadening. These two 
observations are mutually inconsistent. Even though the 
extraction of the in-medium width depends on theory, 
the observed A-dependent reduction of the transparency 
ratio T clearly shows that the co meson is attenuated in 
the target nucleus. This conclusion must be robust. 

Table HVl clearly shows that experimental results have 
not yet converged, and more work is needed to ob- 
tain consistent understanding of the in-medium behav- 
ior of vector mesons. In view of the robustness of the 
method, the in-medium broadening of vector mesons 
deduced from the transparency-ratio measurements are 
hard to rule out. On the other hand, problem(s) have 
been pointed out for all experiments which observed in- 
medium mass shifts, and hence those results must be 
treated with caution and further studies are needed. 
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